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Chapter 1 



Introduction and summary of the 
results 



Today, magnetic fields are present throughout the universe and play an important role 
in a multitude of astrophysical situations. Our Galaxy and many other spiral Galaxies 
are endowed with magnetic fields which are dynamically important [22]. How do these 
cosmic magnetic fields arise? Many astrophysicists believe that galactic magnetic fields 
are generated and maintained by some non-linear dynamo-mechanism, whereby the energy 
associated with the rotation of spiral galaxies is converted into magnetic field energy [Hj. 
The dynamo mechanism only serves as a means of amplification, so that the latter scenario 
requires the presence of seed magnetic fields, whose origin today is still uncertain. 

This diploma thesis is devoted to the study of models and mechanisms which could explain 
instabilities towards the generation of such seed magnetic fields. The structure of this 
report is the following: In chapter 2 we discuss the chiral (abelian) anomaly and related 
questions. The purpose of these developments is to lay the theoretical ground work for the 
calculations in chapter 3. The latter one constitutes the main body of this report and is 
divided into five sections. In sections 3.1, 3.3 and 3.4 we present different mechanisms for 
producing seed fields and derive equations of motion, while sections 3.2 and 3.5 will shed 
some light upon the relations between the individual approaches (see figure [lT]). 
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Figure 1.1: Illustration of the structure of chapter 3, showing the different models discussed and 
how they are interrelated. 

The calculation in section 3.1 assumes that the early universe is filled with a hot plasma 
of charged fermions, whose mass may be neglected, so that chirality is approximately con- 
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served. This plasma is supposed to be in a state of thermal equilibrium and we shall admit 
that there existed a slight asymmetry in the chemical potentials /x/ and /Uj. corresponding 
to fermions of left and right chirality. One may then use a formula - derived by A. Y. 
Alekseev, V. V. Cheianov and J. Frohlich in 1998, [llj - which relates the expectation 
value for the electric current j = ji + jV to the magnetic field B and the difference in 
chemical potentials m — fir- 



By substituting expression (1.1) for the current into Maxwells equations and making some 
simple assumption concerning the charge density, such as = 0, we obtain a system 

of equations which can be solved by means of Fourier transformation. One finds that the 
modes b{t, k) = f d?xe~'^^'^ B{t, x) for 

-> a 

|fc|<^IW-/^r| (1.2) 

grow exponentially with time (a = is the feinstructure constant, we shall set h = c = 1 
througout this report). Hence an equilibrium state with ^ fir is unstable with respect 
to the generation of (electro-)magnetic fields. The reason, why we always talk about the 
generation of magnetic fields is that large electric fields rapidly die out if dissipative 
processes are allowed for [I3] . 

Some of the hypotheses underlying the calculations of section 3.1 seem quite unnatural. 
Chirality is not really conserved, since fermions are massive. The very early universe is not 
really an equilibrium state and the chemical potentials of left- and righthanded fermions 
neither have an unambiguous meaning, nor would they be space- and time- independent. 

In section 3.2 we will show a way out of these difficulties. Based on an analogy with the 
quantum Hall effect we develop a (4+l)-dimentional theory which leads to an effective 
action functional 



= 41 [^^'^] - siSllA^'^] + raA[A(5)|aJ. (1.3) 

In the above formula, ^^^^[A^^)] denotes the (4+l)-dimensional analogue of the Maxwell 
term, ^^^[A''^)] is proportional to the five dimensional Chern-Simons action and the 
boundary term r^Ai^^^^laA] ™ust be introduced in order to assure the gauge invariance of 



the effective action. The equations of motion derived from (1.3) generalise those obtained 
in section 3.1. In particular they yield an equation for the time evolution of the (now 
space-time dependent) difference in chemical potentials. 

While the introduction of a fourth space dimension in section 3.2 could be regarded as a 
mathematical "trick", we take the (4+l)-dimensional point of view more seriously in the 
following section and imagine that there exists a slab A of with L in (4+l)-dimensional 
Minkowski space-time. This slab is filled with massive, (4+l)-dimensional fermions and the 
latter ones are coupled to an external vector potential ^4^^^ Choosing a step function for the 



rr -dependent mass one should then be able to derive the effective action (1.3) up to terms 
O(^), where m is the fermion mass. The massive bulk modes produce the Chern-Simons 
term, whereas FgAi^^^^loA] corresponds to the effective action for the massless, chiral 
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1. Introduction and summary of the results 



boundary modes identified with the (3+l)-dimensional left- and right-handed fermions 
filling the early universe. 

The calculation of det(^^^ -|- im{x^)) for an x'^-dependent mass is rather involved. We 
shall therefore ignore the existence of domain walls in sections 3.3.1 to 3.3.5 and calculate 
the fermion determinant for a constant mass m. To this end one has to compute one- loop 
Feynman diagrams with n = 1, ... ,5 vertices. The vacuum polarisation graph (n = 2) 
yields a contribution which may be absorbed into the Maxwell term Sem by an appropriate 
redefinition of the bare coupling constant (charge renormalization) , whereas the triangle 
graph (n = 3) is shown to produce a Chern-Simons term. The two remaining potentially 
divergent diagrams (n = 4 and n = 5) will not be considered. We end section 3.3 by 
showing how chiral boundary modes occur in the presence of domain walls. 

Equations of motion similar to the ones derived from the (4-|-l)-dimensional theory are ob- 
tained by coupling (3-|-l)-dimensional fermions to an axion field. The axion field approach, 
which appears to be the the most satisfactory, is presented in section 3.4 and yields the 
effective action 

S,ff[A,<p] = SEM[A]-^ JcPiFAF) + J d''x^{d^^<p)id^cP) + W[A], (1.4) 

where (j) denotes the axion field, / some length parameter and Ty[j4] = — ilndet(-P^) is 
the fermionic effective action. 



The action functional (1.4) is equivalent to (1.3) in the case of an x -independent vector 



potential j4^^^ = {Aq, . . . , A4). This can be shown by setting 



Ai = <P, (1.5) 

as will be done in section 3.5. The above identity and the relation £^4 = OqA^ ~ — Hr) 
suggested by the analogy with the quantum Hall effect yield an interpretation for the axion 
field: Its time derivative plays the role of a space-time dependent "difference in chemical 
potentials" between fermions of left and right chirality. 

Another advantage of the axion field approach is that it does no longer rely on any im- 
plausible assumptions, such as "masslessness" of the fermions or a universe in thermal 



equilibrium. One can replace ^^[A] in (1.4) by the effective action for massive fermions. 



Since the equations of motion derived in sections 3.1 and 3.2 may also be obtained from 
the axion field theory, we will restrict our attention to the system of equations involving 
the axion. Of course, the study of these highly non-linear equations is a difficult task 
and the results presented in chapter 4 are only a beginning. We shall try to find special 
solutions and gain some insights into the dynamics by linearising the system of equations 
around these special solutions. 

The form of instability which seems the most appealing to us is the growth of electromag- 
netic fields by parametric resonance. The latter mechanism requires an oscillating axion 
field and is made possible by the presence of a periodic axionic potential f/ Such an 



additional term in the effective action (1.4) is obtained by evaluating the path integral over 
A using a semi- classical expansion based on the stationary phase method. The potential 
U[(f)] might also be useful for finding nontrivial special solutions of finite energy. 
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All the calculations in chapters 3 and 4 are based on the assumption that space-time is flat. 
Physically more relevant results would probably be obtained by considering an expanding 
Friedman-Robertson- Walker universe. While the transcription of our equations of motion 
to the latter model does not appear to present much difficulties, time has not permitted 
us to analyse this new situation and we shall renounce entering into a discussion of this 
subject. 



Chapter 2 

The chiral anomaly 



2.1 Massless fermions and conservation of chirality at the 
classical level 

Before turning to physics, let us present a certain number of results related to the chiral 
(abelian) anomaly - the phenomenon which is at the origin of the mechanisms for producing 
large scale magnetic fields proposed in chapter 3. 

The purpose of this first section is to review some basic notions and to fix the notation. 
We consider fermions in (3+1) dimensions which are described by four-component spinor 
fields 'ip{x). For free particles of mass m the spinor satisfies the Dirac equation 



where the 'j^, /i = 0, . . . , 3 are a set of 4 x 4 matrices satisfying {7^^, 7^^} = 2g^'^ . For the 
metric in (3+l)-dimensional Minkowski space we choose g = diag(l, — 1, — 1, — 1). In the 
chiral representation, the 7-matrices are 



/ = cr'^ is the 2x2 unit matrix and the (T*'s denote the Pauli matrices. The four-component 
field ip may then be written as a bispinor in terms of two-component spinors u and v 



{i-f^'d^ - m)'iP = 0, 





(2.2) 




(2.3) 



and one can introduce the projectors ji and jr such that ^iip = u and 7r'i/' = v: 




(2.4) 



The reason for choosing the subscripts / (=left) and r (=right) will emerge shortly. In 
momentum space the Dirac equation (2.1) reads (j^p^ — m)ip = 0, or more explicitly 



2.1 Massless fermions and conservation of chirality at the classical level 11 



p\ %^)V' = "^V', (2.5) 
where p-a = p^a'^, p-a = p^j^a^ and p = {p^ ,p), p = (p", —p)- If the fermions are massless, 



then (2.5) decouples into two independent equations for u and v: 

{p ■ a)u = f {p ■ (t)u = —pou 



[p ■ a)v = ^ ' I {p ■ a)v = pqv (2-6) 

Since po = \p\ for massless (positive energy) particles and (p-a) is the helicity operator, we 
conclude that u and v describe fermions of left and right chirality respectively. Furthermore 
chirality is conserved, since the dynamics of u and v is decoupled. 

We now couple the massless Dirac fermions to an external electromagnetic field. The 
system obtained thereby is described by the Lagrangian density 

C = -^F^''F^, + tPi^^'{d^,-iA^)il:, (2.7) 
where —^-F^^^Fav is the kinetic energy term for the gauge field and "0 = ij)'^'-^^ denotes the 



conjugate spinor field. In fact the Euler-Lagrange equations of motion derived from (2.7) 
are d^F'^'^ = aip^^ij: and i'^^{d^ — iA^)ip = 0, in which we recognise the Maxwell-Dirac 
electrodynamics. 

According to Noether's theorem, the invariance of a Lagrangian C under global infinitesi- 
mal transformations implies the existence of a conserved current. Given the field transfor- 
mation 



ip{x) — )■ Tp{x) + 6tP{x), 5'4j{x) = f{x)6a 
depending on the infinitesimal parameter 6a, this current is 



and the charge Q = f d?xj^ is then the generator of the transformation (2.8) through the 
Poisson bracket operation 



The Lagrangian density C written in (2.7) is invariant under the local gauge transformation 



V'(x) ^e^"(^)^(x), V(a:) ^ e-^"(^V(a:), A^A + da. (2.11) 



This invariance - for constant a - leads to the conserved current 
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2. The chiral anomaly 



dL 



5a 



■07'^ V'- 



(2.12) 



In addition, the massless theory is invariant under local chiral rotations 



^(^) ^ e^^(^)^V(2;), ^{x) ^ i:{x)e^^^''^^\ A^A + -i^de. (2.13) 



In particular, if is constant the transformations (2.13) are a symmetry of the classical 



system and the corresponding conserved current - called axial current - is 



2.2 Anomalies from the path integral in Euclidean space 
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2.2 Anomalies from the path integral in Euclidean space 

Anomalies occur when the quantum mechanical vacuum functional of a field theory fails to 
have all the symmetries of the classical Lagrangian from which it is derived. The example of 



interest to us concerns local chiral rotations (2.13), which do not leave quantum mechanical 
transition amplitudes invariant. As a consequence the axial current (jg) = {^^^^'^^il)) is 
not conserved for arbitrary external electromagnetic fields. This phenomenon is called the 
chiral anomaly. 

In this section we discuss a technique of anomaly calculation using path integrals in Eu- 
clidean space, the so called Fujikawa method [7] (a useful review on the subject can also 
be found in [8|). As we shall see, the chiral anomaly can be understood as a non-invariance 
of the path integral measure under local chiral transformations on the fermion fields. 

It was mentioned in the previous section that the Lagrangian density C = 4'i'y^{dfj_ — iAn)ip 
leads to the Dirac equation for massless fermions coupled to an external gauge field A. 
The fermion effective action functional S'e//[A] is obtained by performing the integral over 
the fermion fields 

gi5e//[A] ^ I ^^^^^i ! d^xi,i^'^{d,,-iA^)i> ^ ^ = 0, . . . , 3. (2.15) 



S(,ff[A\ is the generating functional for the connected Greens functions of the vector 
currents 

After the Wick rotation 

x° = ^x^ Aq = -iA4, 7° = (2.17) 

the metric becomes g = diag(— 1, —1, —1, —1). This is somewhat unusual and we prefer to 
redefine the 7-matrices according to 

7^^ -i'jf' (2.18) 

in order to obtain g = diag(l, 1, 1, 1). The integral for the effective action in Euclidean 
space then reads 

e-5f„[^] = J p^p^e-i'^'^^^"(^--^^-)'^, (2.19) 
and the Euclidean 7-matrices 7^, fi = 1, . . . , 4 are hermitian, so that 



ipi = ij^'id^ - iA^) (2.20) 



is a hermitian operator with real eigenvalues. Under the local chiral rotations 
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2. The chiral anomaly 



(2.21) 



the exponent in (2.19) is transformed for infinitesimal a{x) as 



(2.22) 



Since 7^ = i7'^7^7^7^ = 7172-^/3^4 anticommutes with 7^, ^ = 1, ... ,4 the axial cmrent 
is invariant under chiral rotations (2.21): j'g = Tp'j'^j^Tp' = ipj^^j^Tp = j^. The rule for 



the transformation of the measure VijjV'ip in terms of the standard Jacobian J reads [T] 



(2.23) 



Note that for Grassmann fields, the Jacobian appears inverted. The latter one can be 
calculated by means of the formula 



detX = e 



trlnX 



(2.24) 



where X denotes some operator. In our case the operators involved are those which perform 



the chiral transformation (2.21 ) on the fermion fields if^ and The Jacobian J is therefore 



J = e 



(2.25) 



If we assume that possesses a discrete spectrum and introduce the eigenfunctions V'n 
corresponding to the eigenvalues An, 



(2.26) 



then the trace in (2.25) can be written as 



tra7^ = ^^('i/'n|a7^|V'n) = J d^xa{x)A{x). 



(2.27) 



The function A{x) appearing in the above formula explicitly reads 



(2.28) 



and is called the anomaly. From (2.22), (2.25) and (2.27) one obtains 



2.2 Anomalies from the path integral in Euclidean space 
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/ 



Vil;Vipe' 



(2.29) 



The right hand side of (2.29) can be expanded to first order in a and the entire equation 
divided by / VipV^pe' ^ x4,pf^ij. 



Using the invariance of the path integral under the change of variable as well as j'^^ = jg , 
this finally yields 



2A. 



(2.31) 



A{x) as it stands in ( 2.28| ) is an ill-defined quantity. We may evaluate it by regularizing 
the large eigenvalues An (which are real, since is hermitian) with a Gaussian cut-off 
and changing the basis vectors to plane waves as 



A(x) = lim '0n(2^)7^^n(3^)e = lim 'tpl^(x)^^e a/2 ip^ix) 

n 



lim , , 

M^oo J (27r)4 



Tr 



327r2 



327r2 



* {F ^F){x). 



(2.32) 



The symbol "A" denotes the exterior product and "*" the Hodge dual. Upon substitution 



of (2.32) into (2.31) we find that in Euclidean space 



^M)A = -.^,*{F^F). 



(2.33) 



The Minkowski-space version of (2.33) is obtained by undoing the Wick rotation: 



167r2 



* {F ^F). 



(2.34) 
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2. The chiral anomaly 



2.3 Chiral currents 

2.3.1 Conserved version of the chiral currents 

The chiral currents jj^ and jr corresponding to fermions of left and right chirality are 
defined as 

jj' = i>j>^^{l - jif = Vi7^^(l + 7')V'. (2.35) 

They are related to the electric current and to the axial current jg by 

f = i^ri^ = + jf, 3^ = ^T^'V' = - if- (2.36) 

The chiral currents and jr are gauge invariant, but not conserved because of the chiral 
anomaly: 

d^jl = T^e^^'^d^A^dpA^. (2.37) 
In certain situations it will be useful to introduce the currents 

~ji,r = fl,r ± ^^e^^'^A^dpA,, (2.38) 
which are conserved, but fail to be gauge invariant. However, the corresponding charges 

4,r- = y d'xf,^, (2.39) 

are not only conserved, but also gauge invariant. More precisely the gauge variation of 
Qz,r amounts to a surface term, which may be dropped. 

2.3.2 Gauge variation of the chiral determinant 



In analogy with (2.15) and (2.16) the functionals T^^^i^] defined through 

^^rA^\ = y p^p^e'^'^'^'^^^''^^'^-^^'')^^^^^')^, (2.40) 
could be regarded as the generating functionals of the Greens functions for the righthanded 



and lefthanded currents respectively. The interpretation of the right hand side of (2.40) 
as det(^^7r,i), as suggested by the formula for Gaussian Berezin integrals [1] 

/ dr\].d% exp[- ^ fjkAk,irii] = detA, (2.41) 

k kl 



2.3 Chiral currents 
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confronts us with a problem though. In fact the determinants of the operators ^^^^,1 
are formally zero. These difficulties may be overcome by redefining e*'"''''^"^' as det{P^^), 
where the operators act on 4-component spinors, but couple the gauge field to posi- 
tive/negative chirality components only: 



Pa = i^{d,-zA,-{l + ^')) 





^ 

^ 

h 



(2.42) 
(2.43) 



In equations (2.42) and (2.43) we used the notation = a ■ d = —a^do — a - V, ^ = a ■ d 
a^do — (7 • V and a subscript A denotes a covariant derivative as usual. 



Under an infinitesimal gauge transformation A 
change as 



A+d9, the generating functionals rr,i[A] 



TrM + de] = FrM] - J d'^()i^)dJJ^- (2.44) 

Hence gauge invariance would require = 0, which is equivalent to the current 

conservation conditions d^{jj^^)A = 0, since 



6rr,i[A] 

dA„, 



(v;7'^-(i±7'Wa = o;Va. 



(2.45) 



But in the preceeding section we have shown that chiral currents are not conserved. A 
theory of massless chiral fermions coupled to an external electromagnetic field is anomalous 
in the sense that it fails to be gauge invariant. 



From (2.36), jr,i = ^{j ± Js) and the conservation of the electric current we find 



"327r^ 



* (F AF), 



(2.46) 



where in the last step we have replaced the anomalous divergence of the axial current 



by its explicit expression obtained in (2.34). The variation of the generating functionals 



rr,/[A] under gauge transformations is obtained by substituting (2.46) into (2.44): 



Tr,i[A + de] = Tr,i[A] T 3^ / ^(^ ^ 



(2.47) 
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2. The chiral anomaly 



2.4 Anomalous commutators 

In this section we would like to determine the equal-time commutators of the current 
components jf.^. The following argument is due to J. Frohlich [13]. It does not claim the 
status of a proof, but hopefully provides a reasonably clear idea about the origin of the 
anomalous commutator and its relation with the anomalous divergence of the chiral cur- 
rents. For a more mathematical approach based on methods of group theory cohomology 
see for example [S] and [TU]. 

Let V denote the the space of configurations of external electromagnetic vector potentials 
A corresponding to static electromagnetic fields. We consider the Hilbert bundle H over 
V whose fibre J-a at a point ^ G V is the Fock space of state vectors of chiral (e.g. 
left-handed) fermions coupled to the vector potential A. 




Figure 2.1: Illustration of the Hilbert bundle Ti, and the gauge transformation U{g) = U{g°'{x)). 

% carries a projective representation U of the group Q of time-independent electromagnetic 
gauge transformations g = {g^{x)) - g^{x) = e*"^^\ a{x) = a{x) independent of x'' = t - 
with the following properties 

1. U{g) : Ta — > ^A+da, 

2. U{g)^lj{x, A)U{g)-'^ = e^'^^^^{x, A + da), 

where ip{x,A) is the Dirac spinor field acting on J^a- The gauge transformation U{g"{.)) 
may be written in terms of the generator G of ^ as 

C/(5"(.)) = e-*"-^. (2.48) 

We used the notation a • G = j d'^xa{x)G{x), where the integration is over space. The 
explicit expression for G reads 



(2.49) 



2.4 Anomalous commutators 
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The first term in (2.49) generates tlie gauge transformation A ^ A + da on A, while 
the zeroth component of the current generates a rotation of the fermion fields. The 
left-handed current appears because we have chosen to consider left-handed fermions. 

Locally the phase factor of the projective representation U oi Q can be made trivial by 



redefining the generators G as G{x) — )• G{x) = —iV 



SA{x) 



+ jf{x,A). 



The operators G{x) generate a representation of the group Q of gauge transformations on 
Ti if and only if 



[G{t,x),Git,y)]=0 (2.50) 
for all times t. We pretend that the right choice for the redefined generator compatible 



with (2.50) is 



G{x) 



-iV 



6A{x] 



(2.51) 



This follows, heuristically, from the fact that jj^ is a conserved current (see section 2.3.1). 
Furthermore, since the current ji^{x,A) is gauge invariant we have 



V-^,jf{x,A) 
dA 



(2.52) 



This enables us to compute the anomalous commutator of the left-handed currents as 
follows (an integration over x and y of the form J dxa{x) J dyf3(y) ... is implicitely under- 
stood): 



6 S 
-iV ■ -J— + jO(x), -iV • -J— + j^{y) 
6A{x) 5A{y) 



0=[G(x),G(y)] = 

= [j?{x),jf{y)]- iV 

= - ^.^'^ii^H^ - y)) i-Mv) + ^^^^■^(|:^(- - y)) i-M^) 



6A{x) Stt^ 



,-^e'^''Ad,Ak{y) - -^e'^''A,d,Ak{x),iV--4 



6A{y) 



[j?{x),jf{y)] + ^^^[B,{x)6{x-y)]. 



(2.53) 



Bi = —e^^^djAk denotes the magnetic field strenght. The calculation for the right-handed 



current is similar. All that changes is the plus sign on the right hand side of (2.51 ) and in 
the subsequent calculations. We therefore obtain the anomalous commutators 



[lr{^),3Uy)] 



i d 
Att"^ dxk 



[Bkix)6{x - y)]. 



(2.54) 



Chapter 3 



The generation of seed magnetic 
fields in the early universe 

3.1 First attempt: Equilibrium statistical mechanics 

This chapter is devoted to the study of several models and mechanisms which could ex- 
plain the generation of magnetic fields in the early universe. In a first attempt we will 
assume that the early universe is a hot plasma of charged fermions and that chirality flips 
constitute a dynamical process slower than the expan sion rate of the universe. Under these 



assumptions the chiral charges Qi^r defined in (2.39) are approximately conserved and we 
may introduce the chemical potentials /i/ and fir canonically conjugate to these conserved 
charges. 

In this first section we shall attempt to show that if there exists an asymmetry in the 
chemical potentials of left- and right-handed fermions, this could lead to the generation 
of large cosmic magnetic fields. 



3.1.1 Current expectation value in thermal equilibrium 

The starting point of our investigations is a formula relating the current expectation value 
in thermal equilibrium to the magnetic field strength, which was obtained by A. Y. Alek- 
seev, V. V. Cheinaov and J. Frohlich in 1998, [11]. We will now go through its derivation, 
that is compute the expectation value of the electric current j in the background electro- 
magnetic field A^. 

The continuity equation for the electric current reads = 0. It can be solved in terms 
of a 3-vector field a: 



f = V-a, j = -doS. (3.1) 

In fact we can define o(x) as 



a(x°, f) = - ^ drjir, x) - j rf'y ^^^ J_ ^ j°(0, y), (3.2) 
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as one may easily verify. A derivation of this result - although not really necessary - is 
given in appendix A. 

The thermal state of the system characterized by the chemical potentials /2 = (/i;, yUr) and 
the inverse temperature /? is given by the density matrix 

^13,(1=^ , T-l-^ - fllQl - HrQr, (3.3) 

where T-L is the Hamiltonian and -Z^ = Tre~^^^ . In this equilibrium state the expectation 
value of the current is 



{j{x))i3,p = -{dQa{x))p^fi = -i{[n, a{x)])i^^ft 



Z, 



-TT{e-^^^ [n^, aix)]} - -^Trje-^^" [wQ, + ^irQr, S{x)]} (3.4) 



The first trace on the right hand side of (3.4) vanishes by cyclicity of the trace. If we 
furthermore use [V ■ a,a\ = 0, that is [j^, a] = 0, we finally obtain 



(W-Air) / d^y{[jUy),^{x)])i3^ft. 



(3.5) 



The commutators of the densities of the left- and right-handed fermions have been calcu- 



lated in (2.54): 



whereas the commutator of the left-handed and right-handed current is zero. Hence 



(3.6) 



my),dkau{x)] = [j°(y) - jT(y), j°(^) + jT(^)] = ^^J^^^i^^^i^ " y)]- (3-7) 



We remove the divergence in (3.7): 



[jl{y),d{x)] = -^B{x)6{x -y)+V^x n(x - y), 



(3.8) 



and substitute this result into (3.5). The second term on the right hand side of (3.8) drops 



out after integration over y. This finally yields 



(i(^))/ 



47r2 



{^il - Hr)B{x). 



(3.9) 
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3.1.2 Equations of motion 



By substituting expression (3.9) for the current expectation value into Maxwells equations, 



we obtain the following system of equations 



V-E = a{f)f3,fl (3.10) 

VxB-doE = a(J)^,^ = -^(w-/"r)^ (3.11) 

V-B = (3.12) 
VxE + doB = 0, (3.13) 

which is supposed to govern the evolution of the electromagnetic field in the early universe. 



The feinstructure constant a = e appears on the right hand side of (3.10) and (3.11) 



because we chose to absorb a factor of e into the definition of the vector potential: 



eA A. (3.14) 

As long as the chemical potentials /U/ and remain constant, the above equations are linear 
in the fields E and B with constant coefficients (at least if we make a simple assumption 
for such as {j^) p^ft = 0). The time evolution can thus be calculated by means of Fourier 
transformation, as will be done in chapter 4. 

We shou ld, h owever, be aware that some of the hypotheses underlying the derivation of 



formula (3.9) appear quite unnatural in the present context. The chiral charges Qi^r are 
not really conserved since fermions are massive. The very early universe is not really an 
equilibrium state and the chemical potentials /i; and /x^ of left- and righthanded fermions 
neither have an unambiguous meaning, nor would they be space- and time- independent. 



In the following sections, we shall try to generalise the system of equations (3.10)-(3.13) 
in order to obtain an equation of motion describing the evolution of a (space-)time de- 
pendent "difference in chemical potentials". This generalisation will ultimately lead to a 
new formulation of the theory in terms of an axion field, which no longer relies on the 
implausible assumptions mentioned above, but could still provide an explanation for the 
generation of cosmic magnetic fields in the early universe. 



3.2 Second attempt: (4+1) -dimensional quantum Hall effect 
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3.2 Second attempt: (4+l)-dimensional quantum Hall effect 



One idea is to imagine that the difference in the chemical potentials fii — fir is generated by 
an electric field E4 pointing in a direction perpendicular to our (3+l)-dimensional world. 
This idea is based on an analogy with the (l+l)-dimensional quantum Hall effect, which 
we shall now discuss. 



3.2.1 The (2+l)-dimensional quantum Hall effect 

It might be useful to summarize some key features of the (2+l)-dimensional quantum Hall 
effect first. This introduction again closely follows the review by J. Frohlich and B. Pedrini 

m- 

A quantum Hall fluid (QHF) is an interacting electron gas confinded to some domain D 
in a two-dimensional plane, subject to a constant magnetic field B^^^ transversal to the 
confinement plane. For D we choose a strip of width L in the 1,2-plane, which is infinitely 
extended along the 1-direction. 

Among the experimental control parameters is the filling factor, defined by 



^(0) 

u = 2n^, (3.15) 

where n^^^ is the (constant) electron density, the component of the mag netic field 
perpendicular to the plane of the fiuid and 2tt the quantum of magnetic flux (in units 
where c = h = 1). 

Transport properties of a QHF in an external electric field are described by the equation 



J(xO,f) =( E{x\x). (3.16) 

In the above formula x is a point in D, j the bulk electric current and E the component 
of the external electric field parallel to the sample plane. Furthermore, aL denotes the 
longitudinal conductivity and an the transverse or Hall conductivity. 

Experimenally, one observes that the longitudinal conductivity, ai, vanishes when the 
filling factor v belongs to certain small intervals. At the same time, the Hall conductivity 
an is a rational multiple of 7^ . Such a QHF is called "incompressible" . 

We will now describe the basic equations describing the electromagnetics of an incom- 
pressible QHF. To this end it is useful to combine the two-dimensional space of the fiuid 
and time to a three-dimensional space-time. The field strength tensor of the system is 
given by 



/ El E2 \ 
Ff,, = -^1 -B \ , (3.17) 
\ -E2 B / 
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where Ei and E2 are the 1- and 2-components of an external electric field and B is the 
component of an external magnetic field, B, perturbing the constant field B^^^ perpendic- 
ular to the sample plane {Btot = B^'^^ + B). We define j*^ to denote the sum of the electron 
charge density at the space-time point x = {x^,x) and the uniform background charge 
density n^^\ From the continuity equation for the electric current density = 
dfj_j^ = 0, the three-dimensional homogeneous Maxwell equations, dF = 0, and from the 
transport equations (3.16) with aL = it follows that 



f = -auB. 



(3.18) 



Equations (3.16), for ul = 0, and (3.18) can be combined to the covariant expression 



{*FY 



(3.19) 



which describes the response of an incompressible QHF to an external electromagnetic 
field perturbing the constant magnetic field B^^\ In (3.19), A denotes the vector potential 
of this external eletromagnetic field {F^i, 



duA^) 



The finite extension of the sample, confined to a space-time region A = D x M, is taken 
into account by setting (Th[ ) to zero outside A, 



(3.20) 



In the above equation, an is the (constant) value of the Hall conductivity inside the sample 
and the characteristic function of the space-time domain A. Taking the divergence of 



(3.19) we get 



e'^'^PdM^F, 



up- 



(3.21) 



Thus d^j^ fails to vanish on the boundary dA of the sample, unless Ff^i,\g\ = 0. For 
arbitrary external electromagnetic fields, there must exist an electric current density jQ\ 
localized on the boundary dA = {(x|x^ = L)} U {(x|x^ = 0)} of the sample space-time 
such that the total electric current density 



Jtot 



(3.22) 



satisfies the continuity equation. These edge currents are chiral, a property which is cor- 
rectly predicted by the naive classical picture of electrons bouncing off the domain walls. 
The current density Ii - localized on the upper boundary = L - is produced by left- 
moving modes (for an appropriate choice of the direction of Btot) and - localized on 
the lower boundary = - by right movers. If If^ denote the corresponding quantum 
mechanical current operators, then the edge currents are given by the quantum mechanical 
expectation value {Xf^)A- 



The effective action 5e//[A] can be found from equation (3.19) relating the current expec- 
tation value to the external electromagnetic field 
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SSeff[A] _ .^^^^ ^ aHe^^'d^Apix). (3.23) 



The solution to the above equation is 



- 5e//[A] = Scs[A] = ~ j d^'xe^^PA^d^Ap = J A A F. (3.24) 

We denote this action functional by Scs - for Chern-Simons action - since it is proportional 
to the integral over the Chern-Simons 3-form. Since the latter one is not invariant under 
gauge transformations of A that do not vanish on the boundary OA of the sample, the 
effective action must be corrected by a boundary term 



S,ff[A] = -Scs[A]+r9A[A\ 



(3.25) 



The boundary term TgAl^bA] is the effective action of the charged chiral modes propa- 



gating along dA. The explicit form of (3.25) and expressions for the boundary currents 
can be found in appendix B. 



3.2.2 Generalisation to (4+1) dimensions 

Let us consider the (2+l)-dimensional Hall sample of width L extended along the 
1-direction. The total Hall current is related to the potential difference V = {^i — ^r) 
between the upper and lower boundary by the formula = crnifJ-i — fJ-r)- If the electric 
field in the bulk of the Hall sample vanishes, then is the sum of two contributions ll 
and supported by so called edge states localized respectively near the upper and lower 
boundary. These edge currents are chiral: if is produced by leftmoving modes and by 
rightmovers, /i; and being the chemical potentials of their respective reservoirs. 

The (l+l)-dimensional system obtained by considering the boundaries of the (2+1)- 
dimensional Hall sample corresponds to a quantum wire in which the left- and rightmoving 
electrons are coupled to reservoirs with chemical potentials ^u/ and fir respectively. The 
total current through the wire is then given by the formula for the Hall current, which in 
turn may easily be derived by considering the bulk of the (2+l)-dimensional Hall sample. 
In fact, if we admit that the entire potential difference between the two edges is generated 



within the bulk of the sample, then (3.16) yields the formula for the total Hall current: 



/I = / dx^j\x'^) = -an / dx^E2{x^) = anil^i - l^r). (3.26) 
Jo Jo 



It seems plausible that formula (3.9) for the current j, which is also the sum of two 
contributions ji and jV corresponding to left- and righthanded fermions, could be "derived" 
in a similar way. To this end we consider a slab of thickness L in (4+1) dimensions, which 
is extended in the 1,2, 3-directions. The upper and lower surfaces represent two copies of 
our (3+l)-dimensional world. Inspired by the analogy with the (2+l)-dimensional Hall 
sample, we place the lefthanded fermions on the top and the righthanded ones on the 
bottom. The potential difference between the two surfaces, which is generated by the 
4-component of an electric field, will be denoted by fii — fir- 
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light-Eimwa 




itmias 



(2+l)-dimensional QH sample 



Slab in (4+1) dimensions 



Figure 3.1: In order to obtain an equation for the time evolution of /i/ — /^^ one imagines that the 
difference in the chemical potentials is generated by an electric field pointing in a direction per- 
pendicular to our (3+l)-dimensional world and exploits an analogy with the situation encountered 
in a (2+l)-dimensional Quantum Hall sample. 



This analogy is illustrated in figure 3.1 It can thus far be summarized by the following 
equations: 

(l+l)-dimensional quantum wire: Massless fermions coupled to an electromagnetic 

field in (3+1) dimensions: 



= anifJ-i - Air) 



J = jl + jr 



(3.27) 
(3.28) 



The above analogies and equation (3.26) suggest to introduce a current density 
J^{x^,x,x'^) in (4+1) dimensions and to proceed as follows: 



(2+l)-dimensional Hall sample: 



(4+l)-dimensional formulation: 



dx^j^{x^,x\x'^) 







cL 



3=1 dx'^ J {x^ , X , x'^) 
Jo 



1 



an! dx^E2(x°,x^,x^) =i/ dx^EJx^ ,x,x^)B(x'^ ,x) 

Jo 47r^ Jo 



(3.29) 
(3.30) 
(3.31) 



Looking at ( |3.29 ) and (3.31) we realise that the (4+l)-dimensional current density J 
satisfies J^{x^ , x, x^) = -^E4{x^,x,x'^)Bk{x^,x). The covariant version of this equation 
reads 



J" 



32^2 



327r^ 



(*(FAF))^ 



(3.32) 



that is 
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J 



327r^ 



* {F AF). 



(3.33) 



The analogue of the field strength tensor in (4+1) dimensions becomes 








El 


E2 


E3 


\ 




-El 





-i?3 


B2 


Vi 




—E2 







-Bi 


V2 




—E3 


-B2 


Bi 







v 


-Ei 


-Vi 


-V2 


-V3 


/ 



(3.34) 



where Vj = Fji denotes some vector field, for which there exists no experimental evidence. 
We shall nevertheless keep the terms in V throughout the following calculations, because 
they will turn out to be related to the axion field which we will introduce in section 3.4. 



In (3.31 ) it has been assumed that the magnetic field does not depend on x^. This condition 
is satisfied if V x F = 0, as will be shown below. Our "derivation by analogy" therefore 
only works in that case. 



3.2.3 Effective action 

Our aim is to find an effective action functional 5ejj[^] from which the equations of 
motion for the electromagnetic field in (4+1) dimensions can be calculated by means of 
the formula 



SA„, 



0, /i = 0,...,4. 



(3.35) 



Of course we know what these equations should be, namely the inhomogeneous Maxwell 



equations for the current density J written in (3.33). It is again possible to find an action 
functional whose functional derivative with respect to the vector potential A yields the 
latter current. The result is the Chern-Simons action 



Scs[A] = d'xe^'''P^''A^d,ApdxA„ = ^ J AaFaF (3.36) 

The apparent solution to our problem is then readily found by adding a Maxwell term 



Sem 



ALa 



d^xF'^'F, 



flU 



(3.37) 



to the Chern-Simons action Scs- The factor of has been introduced in (3.37) for reasons 
of dimensionality and a = denotes the four-dimensional feinstructure constant. In fact, 
setting 



^eff ~ ^EM — Scs, 



(3.38) 
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and introducing (3.38) into (3.35) yields 



^^F''^' = LaJf", = 0, . . . , 4. (3.39) 

An additional set of equations - corresponding to the homogeneous Maxwell equations - 
follows from 



dF = 0. (3.40) 

The latter condition assures that F can be written as the exterior derivative of some vector 
potential A: F = dA. 



However, the action functional Seff written in (3.38) suffers from a serious deficiency: lack 
of gauge invariance. The electromagnetic part Sem obviously is gauge invariant. But the 
Chern-Simons action Scs transforms under a gauge transformation A ^ A + dO like 



Scs[A + de] = ^ ^ d^xe^^'^p^'^iA^ + d^e)d,{Ap + dpe)dx{A^ + d^e) 

= Scs[A] + 
= Scs[A] + 



87r2 J 


1 d' 
'a 


1 

327r2 


/' 

Ja 


1 




327r2 


1 

JdA 



= ^cs[A] + ^ I^^OF A F, A = {Ao,...,A,). (3.41) 

In the above calculation A denotes the volume of the slab in (4+1) dimensions and dA its 
boundary. We have used dF = in the third and Stokes' theorem in the last step. 

One remarks the appearance of a surface term. In order to restore the gauge invariance 
of Seff, it is therefore necessary to add a boundary action Tq\{A\qx) which is not gauge 
invariant either but transforms in a way as to cancel the boundary term coming from the 
Chern-Simons action. 

Luckily we have already examined some possible candidates earlier on. The functionals 
rr,i[A] = -iln J 2?VPV5e^^'^*^'^'^''(^^-^^'*^(i±T'))'/' = -ilndet{p''J), (3.42) 



which have been introduced in section 2.3.2 were shown in (2.47) to transform under gauge 
transformations as 



rr,i[A + dB] = Tr,i[A]T^ J 0{FAF), A = {Ao,...,A3). (3.43) 



The integral in (3.43) is over the boundary dA of the (4+l)-dimensional slab, which 
consists of two copies of a (3+l)-dimensional space-time located at = L and = 
respectively. Defining the boundary term as 
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(3.44) 



we find that the combination — >S'cs'[^]+rgA[^|9A] is invariant under gauge transformations 
in (4+1) dimensions. 



The effective action - which replaces (3.38) - is then 



Seff[A] = Sem[A] - Scs[A] + T9a[A\9a] 



(3.45) 



where the individual terms explicitly read 



Sem[A] 
Scs[A] 



ALa ./a 



+ AF AF 



967r^ ./A 



r^Ai^laA] = -dn{det(^!4U4=i)det(|):4U4=o)}. 



(3.46) 

(3.47) 
(3.48) 



3.2.4 Equations of motion 

The (4+l)-dimensional inhomogeneous Maxwell equations are obtained by substituting 



the action functional Sg// defined in (3.45) into (3.35) and the homogeneous ones are 
dF = 0. The full set of equations then becomes 



5F 
dF 



La 
32^ 
0. 



{F AF)- La 



^Tqa 

5A 



(3.49) 
(3.50) 



Explicitly, equations (3.49) and (3.50) read: 



• Inhomogeneous Maxwell equations, 5F = 

La 



La 
32^ 



*{F AF)- La 



5A 



V • F + 04^4 



-—,V ■ B - La—-- 
An^ 5Ao 



V xB-doE-diV= ^(EiB + VxE)- La^^ 

4vr^ SA 

doE4 - V • y = ^E ■ B 
47r^ 

Homogeneous Maxwell equations, dF = 0: 

V-B = 
V X E + doB = 

VF4 = d^E + doV 
diB = V xV 



(3.51) 
(3.52) 
(3.53) 



(3.54) 
(3.55) 
(3.56) 
(3.57) 



Remarks: 
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In de rivin g the (4+ l)-di mensional equations we assumed that B is x^-independent 
(see (3.31)). From (3.57) we now find that this condition is satisfied if V x 1/ = 0, 



which leads us to postulate that the vector field V can be written as the gradient of 
some scalar field 6: 



V = V(/>. 



(3.58) 



2. For an arbitrary A{x), the boundary action F^a produces a current density on the 
(3+1) dimensional boundary of the slab. However, if A{x) can be chosen such that 



A{x)\9A = 0, (3.59) 

then Fqa does not contribute to the current density. Unfortunately, if B is x'^- 
independent (as it is the case for V = V^), then any solution with non- vanishing 
-B-field is incompatible with the above boundary condition. So we either have to 
relax the condition on V or take care of the currents produced by F^a- 



3.2.5 Projection onto (3+1) dimensions 

In order to obtain the physical quantities observable in (3+1) dimensions as well as their 
equations of motion, one has to project the (4+l)-dimensional fields onto the the (3+1)- 
dimensional space-time identified with the planes x'* = and x^ = L. For arbitrary 
x^-dependent fields it is not possible to express the equations of motion in terms of their 
averaged, x'^-independent counterparts. Some simplifying assumptions concerning the x^- 
dependence of the (4+l)-dimensional fields are therefore inevitable. 

We should also insist on the requirement that the left- and righthanded fermions propagat- 
ing along the surfaces x^ = L and x^ = respectively couple to the same electromagnetic 
vector potential, that is 



^(^)(x,L) = ^(^)(x,0). 



(3.60) 



where A^^^ = {Aq, . . . , ^3) and x = (x*', . . . , x^). This requirement is met if we assume that 
j4(^) = (^(^^,^4) is independent of x^. In this particular case an averaging procedure for 
the fields B, E and V is not necessary. The boundary action for x'^-independent fields is 

Fg^[A(^)] = — ■jln{det(^'i(4))det(^'.(4))} and the boundary current now becomes —^'^^ 



(4) 



-a- 



Furthermore, all the partial derivatives 84, appearing in (3.51)-(3.57) can be replaced by 



zero. The equations of motion in (3+1) dimensions therefore become 
• "Inhomogeneous Maxwell equations" : 



V -E 

V X B-doE 
doEi - V • y 



47r2 5Aq 



(4) 



^(E,B + VxE)-a^ 
4vr2 SA 



La 
4^ 



E-B 



(3.61) 

(3.62) 
(3.63) 
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"Homogeneous Maxwell equations" : 

V-S = 

V X E + doB = 

VEi = d^V 
V X y = 0. 



(3.64) 
(3.65) 
(3.66) 
(3.67) 



Remarks: 



1. The last two equations (3.66) and (3.67) can be satisfied by setting 



V = V(j), E4 = do(j), (3.68) 
where (p is an arbitrary scalar field whose interpretation will be clarified in section 3.4. 



Substituting (3.68) into (3.61 )-( 3.65 ), we obtain the following system of equations 



Inhomogeneous Maxwell equations and equation of motion for (j): 



V -E 

V xB-doE 

ad 



La ^ ^ 5T 
-—.V^-B-a— 



(4) 
OA 



La 
4^ 
La 
4^ 



dT 



(4) 



+ V0 X E)-a—^ 
6A 



E-B 



Homogeneous Maxwell equations: 
V-S = 

V X ^ + = 0. 



(3.69) 

(3.70) 
(3.71) 



(3.72) 
(3.73) 



2. At the beginning of this chapter, we set out to derive an equation describing the 



time evolution of the difference in chemical potentials appearing in (3.11). We now 
show that the system ( 3.61| )- (|3.67[ ) indeed provides such an equation and therefore 



generalizes the system (3.10)-(3.13) obtained in section 3.1. Setting V = 0, neglecting 



the term in F^'^ and introducing the notation 



IJ'l — fJ-r = —LE4 
consistent with the (4+l)-dimensional QH-analogy we find 

• Inhomogeneous Maxwell equations and equation of motion for fii — Hr 

V-E = 



(3.74) 



V AB-doE 

doifll - fir) 



a 

al? 
'47r2 



(/i/ - Hr)B 

E-B 



(3.75) 
(3.76) 

(3.77) 
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Homogeneous Maxwell equations and additional condition for /i/ — 



V-B = 

VAE + doB = 

V{fil-fir) = 0. 



(3.78) 
(3.79) 
(3.80) 



Note that the condition of x^-independence for E and B may be relaxed if the field V 
is ignored. Indeed, a set of equations similar to (3.75)-( [3^ ) is obtained by defining 

and the dif- 

L 



the electric field in (3+1) dimensions by E{x^, ^) = Jf^ dx'^E(x 



ference in chemical potentials by — /i^ ^ — /g dx'^E4{x'^, x, x^). However, we shall 
not pursue this possibility any further and content ourselves with the phenomena 
produced by x^-independent fields. 



3.3 Massive (4+l)-dimensional fermions confined to a slab 
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3.3 Massive (4+1) -dimensional fermions confined to a slab 



After a first, more intuitive derivation of the (4+l)-dimensional theory based on an analogy 



with the quantum Hall effect, we will show in this section that the effective action (3.45) 



may also be obtained in another way. While the introduction of a fifth dimension in the 
preceeding discussion could be regarded as a mathematical "trick", we now take this 
point of view more seriously and identify the (3+l)-dimensional world with a domain 
wall in a (4+l)-dimensional space-time. More precisely we admit the existence of a slab 
A = (x|0 < < L) filled with (4+l)-dimensional fermions of mass m. Two domain walls 
at = and = L are introduced by setting the fermion mass to —M outside A. The 
calculation of the fermion determinant in the limit M — )• oo produces the effective action 



(3.45) up to terms of order O(^) and the massless chiral fermions in (3+1) dimensions 



appear as surface modes localized near the domain walls. 



3.3.1 Calculation of the fermion determinant 

For the moment we ignore the existence of domain walls. Our aim is to calculate the 
fermionic effective action Seff[A,m], which is defined as 



Seff[A,m] = —i\n / Vif^V^l^e 



J J d^x{tp{iPA-m)ip} 



-ilndet(^^ + im) 



(3.81) 



In the present context the symbol iPa denotes the Dirac operator in (4+1) dimensions, 
iPy^ = ij^idfj, — lA^), /X = 0, . . . , 4, and 7^ = i'j^, so that {7^, 7"^} = — 2 = 2g^^. In order 



to evaluate the fermion determinant (3.81) we define the operators^ and 5 through 



/|x) 
{x\S\y) 



{x)\x)=j^A^{x)\x) 
d^k ^ + m 



(27r)5 A;2 - tv? 



-ik(x—y) 



(3.82) 
(3.83) 



Since {x\{ip — m)S\y) = {ip^ — m){x\S\y) = 5{x — y) = {x\y), the following identities 
hold 



i{P-m)S = 1 
iP + im){l + 5/) = Pa + im. 



(3.84) 
(3.85) 



Plugging ( |3.85| ) into ( |3.8ip yields Seff[A,m] = -i\n[dei{P + im)det(l + 5/)]. The first 
determinant is independent of A and adds an irrelevant constant to the effective action, 
which we shall ignore. Using det(l + S^) = expTrln(l + S^) and writing the logarithm as 
a power series we find 




(3.86) 
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The explicit expression for Tr{S^)"' reads 



Tr(S'/)" = tr y" dx{x\S4 . . . S4\x) 

= tr y dxi . . . dxnS{xi - X2)4ix2) ■ . .S{xn - xi)/(xi) 
dxi . . . dxnA^-^ (xi) . . . j4^„ (x„) J ^2^-j5 • • • "^271^ 



X tr 



pi 



k"^ — m? 



-ik-i{x\—X2) 



(3.87) 



where tr denotes the trace over the 7-matrices. An ie-prescription is implicitely contained 
in the definition of the fermion mass. We shall therefore write —m?' rather than — + ie 
in expressions such as (3.87). A notation consistent with the Feynman rules could be 
obtained by writing = Spi—iA)- The right hand side of (3.87) is then identified as the 
sum of the contributions of one-loop diagrams made of n Feynman propagators Sp = iS 
and n vertices —i^. 

The contributions from the diagrams with n < 5 are potentially divergent, while those 
from higher order diagrams vanish in the limit m — >■ cx). We will now proceed to the 
calculation of the amplitudes for n = 1, 2 and 3. 



3.3.2 Contribution in A 



Figure 3.2: Diagram corresponding to w^. 



Tr(S(A) is easily seen to vanish. We will nonetheless treat this simple case systematically 



in order to warm up for the following calculations. From (3.87) we find for n=l 



TV(54) = j d'xA,{x) I ^tr 



+ m 



k^ 



(3.88) 



The A;-integral gives the amplitude corresponding to the diagram in figure 3.2 



d^k 



rtr 



^ + m 
J^2 _ 



(3.89) 



We now look at the different contributions found by evaluating the numerator in (3.89). 



The term linear in m vanishes since tr7^ = and the remaining one is is odd in k, so 
J d^kti[...] yields zero. Hence Tr(54) = 0, so that 



q("=i) _ n 



(3.90) 
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3.3.3 Contribution in 




Figure 3.3: Diagram corresponding to lu 



The calculation of the amplitude associated with the vacuum polarization diagram shown 
in figure 3.3 is somewhat more involved. In the limit m — ?■ oo it will lead to an expression 
of the form J d^xF'^, with a divergent prefactor. This contribution may be combined with 
the Maxwell term to yield a finite result, if the bare charge of the fermions is appropriately 
redefined (charge renormalization) . 



For the actual calculation we proceed as before. Starting from equation (3.87) for n 
and changing variables according to ki ^ k, k2 ^ k — p one finds 



TV(54)2 = dx dyA^{x)A,{y) 



d^p f d^k 



(27r)5 J (27r)5 
The next step is to define the function uj^'^ as 



tr 



7 



-ip{x-y) 



(3.91) 



d^k 
(2^ 



rtr 



+ m 



A;2 



— j( + m 



171"^^ {k — pY 



(3.92) 



which corresponds to the amplitude associated with the vacuum polarisation diagram 
shown in figure 3.3 The integral (3.92) defining uj^^ - which looks quadratically divergent 



for large internal momentum k - may be given a meaning using Pauli- Villars regularization. 
This amounts to coupling the gauge field to additional spinor fields with a very large mass 
XsTTi and (eventually) different statistics. Such a prescription implies the replacement 



UJ' 



{p,m) 



{p,m)+^CsUj^"'{k,Xsm), 



(3.93) 



where the substitution is understood under the integral sign in (3.91). The constants C. 



will be chosen in order to remove the divergence of the integral. The following calculation 
is based on the discussion of the analogous problem in (3+1) dimensions found in the book 
by Itzykson and Zuber [IJ . 



Denoting the large masses Xsin collectively by A we find 



^Chapter 7: Radiative corrections, p. 319-323 
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u;'^'^(p,m,A) 



(M 



rtr 



^ + m ^ — ^ + m 



k'^ — rrfi^ {k — pY — rrfi 



+ ^ Cs{m — > Asm) 



s=l 



d^k { kt'ik - pY + k^ik - pY - g>"'{k'^ -kp- m^) 



(2^)^ 



(/c^ — m?)[{k — pY — rn?) 



S 1 

+ ^Cs(m Asm) \ 

s=l J 



(3.94) 



By rewriting denominator of (3.94) using the formula [T] 



^2 _ 



Jo 



(3.95) 



and introducing auxiliary five- vectors zi and Z2, one can generate the integrand in (3.94) 
by differentiation 



ujf,y{p,m,K) = j dai J da2 J 



d d d d 
+ 



dz^ dz2 dz^ dz2 



d d 



^ gi(ai(fc2_m2)+Q2((fc-p)2-m2)+2ifc+22{fc-p)) _|_ ^(^^(^ _^ X^rn) > 

S = l J 21 = 



=22 =0 

(3.96) 



Integrating over k (using Fresnel) and performing the required derivatives yields 



uj^y{p,m,K) 



ITT POO fOO 



g 4-" daida2 
^7ri Jo Jo (ai + 02)^ 



2aia2P^Pu 
{ai + a2Y 



aia2P 



+ m" 



t|— m^(ai+a2)+ "^"^ 



+ ^Cs(m Asm) y 
s=l J 



(ai + a2Y oi + "2 

(3.97) 



The polynomial in p appearing in the above integral can be rearranged to read 



2\ OL\a2 



{a\ + a2Y 



2 «i«2 2 
m — —p 



a\ + 02)^ a\ + 02 



(3.98) 



We will use the symbol Aw to denote the contibution to uj^u{p^m, K) coming from the 
second term on the right hand side of (3.98) and the analogous terms involving Asm. Since 



Pauli-Villars regularization preserves gauge invariance and the latter term does not seem 
to exhibit this property, we shall admit that its contribution is zero. Under the assumption 
that Acj vanishes, uj^y{p,m, K) reduces to 
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u}f„y{p,m,A) = —{p,j.Pu - g^uP^)uj{p^,m,A), (3.99) 

An 2 

where uj{p'^,m, A) is defined as 

^(p2,m,A) = r rdaida2^^^f:C.eH-"^(°^+"^)+^^1. (3.100) 
■^0 Jo (ai + 02)2 

We used the notation Co = 1, rriQ = m? and = A^m^. Introducing a factor 1 = 
dp6{p — ai — 02) under the integral sign and changing the variables according to 
Ui = pf3i we find 



uj{p'^,m, A) 



JO JO p2 — 



^ ^^gip(-m2+/3ifep2) _ (3.101) 



Since < ^i,/32 < 1 and /3i + /32 = 1 it follows that /3i/32 < j. So if we pick < 4^2^ 
then — /3i/32P^ is positive and the integration contour in the complex p-plane can be 



rotated by — | in such a way that the integral !„ defined in (3.101) reads 



s 

Ip = lim Cs 

s=0 



2i 

r 
r/2 



(e*4 - 1) + 



(-0 



dp 
p2 



(3.102) 



see figure 3.4 The first term on the right hand side of 3.102 is the contribution coming 
from the small quarter circle of radius r] (up to a correction of order 0(y^)) and the 
second term corresponds to the integration over the negative Ini/o-axis. 



Figure 3.4-' Rotation of the integration contour in the complex p-plane. 



Changing the integration variable from p to p\/m1 — /3i/32P^ and integrating by parts, we 
obtain 
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S r ^ . I n, 

"^"^ Lr/2 (-i)5 Jr, 



dp 



S 

lim V Cs 

s=0 
S 

lim V Cs 

s=0 



2i 



- 1) + -^^ml-^,P2vH\ 



(-0 



dp 



p2 



2i 
r 

7^2 



-z 2 



(3.103) 



The constants Cs have to be chosen such that the divergent term proportional to 2 
disappears. This is achieved by setting 



5 S 
s=0 s=l 



(3.104) 



Introducing the remaining term into (3.101) then leads to 



a;(p^, m, A) 



^y2Cs [' dPPil - /3)Vm2-/3(l-/3)p2. (3.105) 

-i)2 „_n ^0 



Developing (3.105) in powers of p yields the result 



uj{p'^, m, A) 



3{-iy. 



S 

m 



(3.106) 



where the O(^) also depends on the constants Cg and A<j, but unlike the first term - which 
corresponds to 6t;(0,m, A) - does not diverge if the regulator masses are taken to infinity. 
We now introduce the ultraviolet cutoff A through the relation 



S S 
s=0 s=l 



(3.107) 



Given the value of A, a possible choice compatible with (3.104) and (3.107) would be 
5 = 2, Ci = -2, C2 = 1, Ai = A and A2 = (A - 1) + 2A. 



Upon substitution of (3.107) and (3.106) into (3.99) we find 



iA\m\ 



+ .... 



(3.108) 



The additional contribution coming from the O(^) in (3.106) has been indicated with 
triple dots since it vanishes in the limit m — )• cxd. By introducing the leading term into 



(3.91 ) one obtains 
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dx / dyA^{x)At,{y) 



-ip{x-y) 



dx / dyA^{x)A^{y) I 2g- 



d d d d 



iA\m\ 
' 247r2 
iA|m| 
' 247r2 

/ d'y I d'x{d^A''{x)-d''A^^{x))id^A,{y)-d,A,iy))6ix-y) 

iA\m\ 
' 247r2 



5 5 



(9ya dxn dyu dxy dy. 



- y) 



(3.109) 



The contribution to the action Sef / coming from Tr(S/)2 is fomid from (|3.109|) and (|3.86|). 
It reads 



.(n=2) 
'e// 



A|m| 
487r2 



(3.110) 



As we shall see in section 3.3.5, it is possible to redefine the bare charge appearing in 
the Lagrangian density in such a way, that the S^^j'^'^ and the Maxwell term with bare 
coupling constant {S^em ~ ~4Zc? / d^xF"^) add up to Sem = ~iZa f d^^^"^- 

3.3.4 Contribution in 




Figure 3.5: Diagram corresponding to 



Tr(S(4)^ leads to the Chern-Simons 5-form A A F A F. Equation (3.87) for n = 3 and after 
a change of variables k2 ^ k, ki ^ k + q, ^ k — p reads 



d'x J dS J d'^zA^{x)A,{y)A,{z) J 
d^k 



d^p f d^q 



(2vr) 



rtr 



It + q + m 
{k + g)2 — im? 



r 



(27r)5 J (27r)5 



^2 _ (J^ _ p^2 _ j^2 



(3.111) 



The function uj^^^P{p, q) which (up to a factor of —1) corresponds to the amplitude associ- 
ated with the triangle diagram shown in figure |3.5| is 



a;^-^(p,(?) 



d^k 



rtr 



{k + qy 



fe2 — m"^^ {k — pY 



(3.112) 
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Next, one tries to identify those terms in the nominator which do contribute to the final 
result. The term in vanishes since tr(7'^7'^7'^) = 0. Those in m are necessarily of the 
form tri^-ft'Yh^) = 4.e°""'^^aab,3 (7^ = 17^) since the nonvanishing traces of this type 
must contain exactly one matrix 7^. Furthermore the contributions linear in k cancel and 
we cannot have a = b = k because the e-tensor is completely antisymmetric. The only 
remaining possibility is a = q and b = —p that is 

tr(^7'^7"(-^)7') = -4e"^'''''qaPp. (3.113) 

The sum of the terms in m? and does not seem to meet the requirements of gauge 
invariance and Lorentz covariance, so we shall simply admit that this contribution vanishes. 
Under this assumption we find 



u^^'^P, q) = -4me"^^^^g„p;3w(p, q, m), (3.114) 
where the function uj{p, q, m) is defined as 



uj{p,q,m) 



d^k 1 

(2^)5 ((A; + g)2 - m2)(A;2 - m'^){{k - pf - m?) 
^1 dxl dyj (2,r)5 [k^ + 2kP + M^r ^'''''^ 



JO 



The variables P = qy — px and M'^ = cp'y + p'^x — m? have been introduced for simplicity. 
The integral over k is convergent and could be performed using one of the formulas listed 
in the appendix of p|. It might nonetheless be useful to do this calculation in more detail, 
were it only for earning a better understanding where such formulas come from. 



Changing variables in (3.115) according to A; — )• k + P (which is legitimate for a convergent 
integral) and remembering the ie-prescription implicitely contained in the definition of m 
we obtain 



/ 



d^k 1 f d^k 



(27r)5 [fe2 + 2A;P + M2]3 J (^2Tr)^ [k^ - {q^ {y^ - y) + p'^ (x^ - x) - 2qpxy + m'^ - ie}f' 

(3.116) 



Since k = kf. — k\ the poles in the complex fco-plane are located at 



ko = ±[P + {q^iy"^ - y) +p^{x^ - x) - 2qpxy + m^ - ze}]5. (3.117) 



They are illustrated in figure [3^ by black dots. The integration contour can thus be rotated 
in the complex A:o-plane as indicated in the figure: 
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Figure 3. 6: Rotation of the integration contour in the complex fco-plane. 



27r J (27r)4[fcQ2_^2_{. . .}]3 

d^k' 1 

oo . (2vr)4 



~ dk'^ 



2-K 



1 



(27r)5 [fe'2 + {...}]S 



(3.118) 



where in the second step we made the change of variables k'^ = iko, k' = k (Wick rotation) 
and {. . .} stands for the expression in between the curly brackets appearing in (3.117). 
This rotation being performed, we may again ignore the ie-term. 

The last expression is an integral in 5-dimensional Euclidean space which may be evaluated 
using the formula [J| 



1 _ 1 b(t-")r(a-f) 
(27r)« {k^ + 62)a - f(a) ' 

We therefore obtain a relatively simple expression for the integral over k 
d^k 



1 



1 



(2^)5 [fc2 _{...}]3 647r2 y/g2(y2 _ y) + ^2(^2 _ _ 2qpxy + ^ 



(3.119) 



(3.120) 



and from (3.120), (3.116) and (3.115) 



Lo{p,q,m) 



327r2 



1 pl—x 

dx / 

Jo 



\/ 9^(2/^ — y) + p'^ix'^ — x) — 2qpxy + 



^ 1 ^/ 1 N 

+ o(-^). 



647r2 



(3.121) 



Substituting the above expression for u into (3.114) yields 
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oo^''Hp,q) = -j^^^e'^^'^'qaPp + O(^), (3.122) 
Idvt^ \m\ 



thus in the limit m — )• cxd we obtain 



Tr(5'4)3 = _^sgn(m) j d^x j d^y j d'> zA^{x)A,{y)As{z) 



(27r)5 J (27r)5 

^sgn(m) j d'x J d% J d'ze'^''^^'A^{x)A,{y)As{z)^^5{x - y)6{x - z) 



I 



167r2 

—i 



sgn(m) j d^xe'^'^-'^^A^daA^dpAs 



„sgn(m) / AaFAF. (3.123) 
647r^ 



The contribution of Tr(5yi)'^ to the effective action S^fj is found from (3.86). It equals 



half the Chern-Simons action defined (for positive m) in (3.36) 



Si^f^ = -lsgn{m)Scs = -r^^M^) I ^AFAF + 0(-^). (3.124) 



3.3.5 Charge renormalization 

Infinities like the one encountered in the calculation of the vacuum polarization graph 
teach us that the parameters such as charge, mass, . . . appearing in the Lagrangian are 
not necessarily observable quantities. The Lagrangian is written in terms of the bare 
quantities tp^, and e*. In our case 

£ = i^'ii^id^ - iA,) - m'}i;' - J-^F^,F'^^ (3.125) 

where = (e^)'^ is the bare feinstructure constant. The reason why we do not introduce 
a bare gauge field A'' is that the condition of gauge invariance requires the combination 
— iAfj^ to appear in the Lagrangian. 

If the theory is renormalizable, then the bare quantities can be related to the physical 
ones by renormalization constants ^, /u and Z as 

^|;' = ^'2{A)^, m'' = MA)m, = (3-126) 

where the renormalization constants are functions of the ultraviolet cutoff A (for a short 
treatment of this subject see for example the book by Collins-Martin-Squires |^ [3]). The 



Lagrangian (3.125) may then be expressed in terms of the physical quantities as 
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C = i{K)i^{i^^'{^^ - lA^,) - ^l{K)m]i; - Z(A) 



4Lq 



(3.127) 



Each of these renormalization constants can be written as a series in the effective couphng 

,2 _ J_ 
"""^ 137' 



a 



e(A) = 1 + 6(A)« + 6(A)«'.-- 
/x(A) = l + /ii(A)a + /i2(A)a^.. 
Z(A) = l + Zi(A)a + Z2(A)a2 + ... 



(3.128) 
(3.129) 
(3.130) 



The functions /i„ and Z„ contain the divergences for A — )• oo, which are thus absorbed 
into the definition of the bare quantities in such a way that the physical quantities remain 
finite. The theory is cahed renormahzable if this can be done while keeping the form of 
the Lagrangian the same as the original. 

As a consequence of these remarks, the fermion mass m appearing in the calculation of 
the vacuum polarisation graph should be replaced by ^(A)m and the propagator S by 
^(A)5. However, the most divergent term corresponds to order zero in the power series 
developments of /z and ^, (that is /i(A) ~ 1 and ^(A) ~ 1) and therefore reads 



g{n=2) 



(3.131) 



as in (3.110). This term coming from the calculation of the vacuum polarization graph 



may be cancelled by an appropriate choice of the coefficient Zi{A), which determines the 
first order correction to the bare coupling constant a''. In fact, by adding the last term in 



(3.127) to (3.131) and developing Z{A) to first order in a, we obtain 



4La 



(1 + Zi(A)a + . . .) J d^xF^^F^^r (3.132) 



Since the above sum should produce the result Sem — — JEa 

/ (fxF^^,F^"' this yields the 

identity 



Alml 



1 



487r2 4La 
which determines the the coefficient Zi(A): 



(l + Zi(A)a + ...) 



1 

'iZ' 



(3.133) 



Zi(A) 



LK\m\ 



(3.134) 



The bare coupling constant a'^ corrected to first order in the physical coupling a therefore 
reads 



a" = a{l + 



LA|m| 
T2^ 



-a + 



(3.135) 
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3.3.6 Massless boundary modes 

If the massive (4+l)-dimensional fermions are confined to a slab of thickness L, then 
besides the contributions from the heavy bulk modes - which yield the Chern-Simons 
action Scs - there exist massless chiral boundary-modes localized near the edges at = L 
and = 0. These edge states are identified with the left- and righthanded fermions in 
(3-|-l)-dimensions. The effective action of these chiral modes corresponds to the term Tq\ 



appearing in (3.45). 



The confinement of the (4-|-l)-dimensional fermions to the slab could be modeled by 
choosing the fermion mass m{x^) as 

f 4\ f m , < < L 1QRN 

"^(^ ) = I -M , -L^<x'<0 ' (^-^^^^ 

and identifying the (3-|-l)-dimensional surface x^ = —L^ with x^ = L (periodic boundary 
conditions). M is eventually taken to infinity. The (4-|-l)-dimensional space-time is thus 
divided into two domains, 



ni = {x,-L^ <x^ <0), 



^2 = (x,0 < < L), 



(3.137) 



I I f5) A 

see figure 3.7 Unfortunately, the calculation of the determinant det(^^(5j + im{x )) is a 
difficult task. Papers related to this problem are [14j and [15j . but they mainly discuss the 
question of anomaly cancellation. An explicit calculation of the effective action in (2-1-1) 
dimensions for a single domain wall can be found in |16j . 

We shall content ourselves with the calculations of sections 3.2 - 3.4 for an x'^-independent 
mass, which are valid within the domains ili and as long as one does not approach the 



domain wall too closely. Making the simplifying assumption that formula (3.124) remains 



valid up to the domain wall, we may again deduce the form of the boundary action Tq\ 
by looking at the gauge variation of the Chern-Simons term. Note that the sign of S'^J^^ 
depends on the sign of the fermion mass. For m(x^) defined in (3.136) we therefore obtain 



1 oidni) 



[A] 



1 nidn2) 

2^cs 



(3.138) 



where 

transformation A ^ A + d6 is then found from (|3.41|): 
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1,2. The variation of si"f^'^ under a gauge 



'eff 



si;f\A + de] = sl]j'>[A] + 



(n=3)r 



s: 



(n=3) 



1 

32^ 
1 

32^ 



1 

e{F A F) 



e{FAF) 
1 



327r^ 



9{FAF). (3.139) 



Looking at (2.47), we realise that gauge invariance may be restored by adding the boundary 
action 
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Tan, [Alan,] = r,[^U4=o] + TiIAI^^l] = -ilndet{(P^U4=i)(P^U4=o)}. (3.140) 



This is again the sum of the effective actions for lefthanded (3+l)-dimensional fermions 
locahzed at = L (r;[^|3,4=x,]) and the effective action for righthanded fermions locahzed 



at = (rr[^|^.4^Q]). The action functional (3.45) - multiphed by a factor of ^ - and 



the corresponding equations of motion (3.61 )-( 3.67) may be obtained by an appropriate 



choice of the length parameter I appearing in the Maxwell term —^[^ f F"^- Indeed, let 
us renormalize the charge in such a way, that the Maxwell term and the contribution of 
the vacuum polarization diagram yields f inside ^2- The coefficient multiplying 
the corresponding term in $7i will then be much larger, since M ^ m. This is equivalent 
to saying that the Chern-Simons and boundary currents associated with ^2 (both inde- 
pendent of the fermion mass) will be strongly suppressed. Their contribution may thus be 
neglected when calculating the average over x^. Since furthermore ^ ^ L, we obtain - 
after averaging - the equations of motion (3.61 )-( 3.67) if / = 2L. 



Figure 3. 7: Step function for the x'^-depen- 
dent mass and definition of the domains fJi 
and ^2- 



Figure 3.8: Massless chiral zero mode loca- 
lized near the domain wall at = 0. 



The eigenfunctions of the 5-dimensional Dirac operator corresponding to the chiral boun- 
dary modes can be calculated (in the limit m, M — )• oo) if the 5-component vector potential 
^4^^) is independent of the coordinate x^. In this case the five dimensional Dirac operator 
can be written as 



- m(x^) = IPa - l^di + il^A^ - m(x^), (3.141) 

where A = (Aq, . . . ,^3) and denotes the four dimensional Dirac operator. We are 
looking for solutions ^'^ to the Dirac equation {ip^'^^ —m{x ))'^r = 0, which are localized 
near the domain wall at = 0. To this end we write = '^rix)r]{x, x'^), where x = 
{x^, . . . ,x^) and ipr describes a righthanded fermion in (3+1) dimensions, coupled to the 
exterior electromagnetic vector potential A: 



iP^lPr = 0, J^iJr = tpv (3.142) 

The resulting equation for r] is not easy to solve (unless A^ = 0) and we shall content 
ourselves with an approximate solution, which becomes exact in the limit m, M — )• co. We 
propose the ansatz 7/(x,x'^) = N{m, M)e~ ■l^o dym{y)-iA4x ^^j. ^^iq function t], so that 
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^r{x,X^) = lPr{x)N{m, M)e- ^0 dym{y)-iA^x^ 



(3.143) 



N{m, M) is a normalization constant. Applying the five dimensional Dirac operator to the 
function ( |3l43l ) we find, using ( [3l42| ) and ( |3l4T| ) 



[iP 



(5) 

A(5)(z) 



(3.144) 



4 

The function |t/P = A^(m, M)^e~^-^o ig localized near the domain wall at = 

and tends to the delta-function 5(x^) in the limit m, M — )• oo. Hence = 'i/'r?? is localized 
near x^ = and satisfies the Dirac equation in the latter limit, as x'^ri{x,x'^) — )• (see 



figure ?? for an illustration) 



The effective action for these righthanded boundary modes (in the limit m, M — )• oo) is 
formally given by 



(3.145) 



On the right hand side we recognise the effective action (2.40) for righthanded (3+1)- 
dimensional fermions. 
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3.4 Third attempt: Axion electrodynamics 
3.4.1 Effective action 

A similar set of equations as the one derived from the (4+l)-dimensional theory, is obtained 
by coupling the (3+l)-dimensional fermions to an axion fielcQ As it will turn out, the time 
derivative of the axion field corresponds to what we called fj-i — fir in section 3.2. 

For a short review on the subject of axions, see for example [17J. We will consider here 
the so called model independent axion, first described by Witten [18J. In string theory 
there appears an antisymmetric tensor field -B^jy, which in (3+1) dimensions possesses one 
physical degree. The associated field strength is not just the curl of B, but is made gauge 
invariant by adding a Chern-Simons 3-form 

H = dB-^{AAF), (3.146) 

and therefore 



dH 



327r' 



;iFAF). 



(3.147) 



Formula (3.146) is valid in a flat space-time (i? = 0) and for a system in which the 
electromagnetic field is the only gauge field present. We shall accept these results without 
further justification. For more information see the book by Collins-Martin-Squire^ [3]. 

The equation of motion is 



6H = 0, (3.148) 
where 5 is the co-differential. If we write Y = *H, this implies 



*dY = 6H = 0, (3.149) 

or simply dV = 0. The latter equation is solved by setting Y = da. In (3+1) dimensions, 
y is a 1-form, so a is a pseudo-scalai|^ We now define the axion field (p by 



a=^, (3.150) 
la 

where the parameter I - with the dimension of length - has been introduced in order to 
obtain [(/>] = Hence the axion field 4> is related to the dual of H by 



few days before the completion of this report, we reahsed with some deception that the results 
derived in this section are not really new. More than a decade ago, M. S. Turner and L. M. Widrow have 
proposed that the axion field could provide a source term for large-scale magnetic fields, see |22| . 
^Chapter 10: Str ing Th eories, in particular section 10.7 on Anomalies. 



^From equation (3.1521 below it follows immediately that is a pseudo-scalar, since □ is a scalar and 



*{F /\ F) ^ E ■ B a, pscudo- scalar. 
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1 

la 



dcf) = *H. 



(3.151) 



We shall now explain why is an axion field. Applying the co-differential 5 to equation 



(3.151) we find 



la6*H, that is □(/> = l*dH. Using (3.147) the latter equation yields 



la 
32^ 



(FAF), 



(3.152) 



which is the Euler-Lagrange equation of motion corresponding to the action functional 



(3.153) 



In the second term we recognise the standard coupling of an axion to the gauge field A. 
This term may be understood as arising from coupling fermions to the axion as 



(3.154) 



A system of charged fermions coupled to an electromagnetic vector potential A and to an 
axion field (j) is therefore described by the action functional 



(3.155) 

The / in front of the dfj_(j)'ipj'^'y^ip term follows from dimensional consideration, whereas the 
factor of ^ must be introduced in order to obtain the equation of motion (3.152 ). Carrying 



out the integral over the fermionic degrees of freedom we find the effective action 



JS^ff[A,4 



J V^Vi^e^'^^'^^^'^^ = det(P^^.,^^^, 



(3.156) 

In connection with the chiral anomaly we have seen in section 2.1 that 



^^KPA+Ld,W^) = e-i^^^(^^^)det(^^) 



(3.157) 



Hence the effective action may be written as 



S,ff[A,cP] = W[A] 



32^2 



0(F A F) + / d^x{ 



4a 



-F'^^'F, 



(3.158) 



which is expression (3.153) up to the fermionic effective action VF[^] = — ilndet(P^) 
and the Maxwell term. Note, that it is not necessary in this approach to assume that the 
fermions are massless. The same calculations go through in the case of fermions of mass 
m, except that the fermionic effective action then reads = — ilndet(P^ + im). We 

shall henceforth consider this more realistic situation. 
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3.4.2 Quantum fluctuations and cixionic potential 

A transition amptitude from a configuration {<j)in,Ain) of the electromagnetic and the 
axion field at some very early time ti to a configuration {(pout, Aout) at a much later time 
t2 can be computed from the Feynman path integral 



(3.159) 



with boundary conditions {(pih) , A{ti)) = {(f)in,Ain) and ((/)(t2), ^(^2)) = i4>out,Aout) (the 
term W[A] in Seff[A, (p] also depends on the boundary conditions at times ti, t2 imposed 
on the fermion fields, which have been integrated out). 

The integral over the gauge field configurations yields the effective action for the axion 
field, 5e// [</>], so that 



(3.160) 
(3.161) 



We evaluate the integral (3.160) by using a semi-classical expansion based on the stationary 



phase method. The equation for the saddle point is 



6S. 



effl 



(3.162) 



and we shall denote the solution of (3.162 ) by (pci- From (3.161 ) one then finds the following 
expression for Seff[(l)d] 



S,ff[(P,i] = -iln J VAe'^^ff^^' 



<f>cl\ 



(3.163) 



Again, we use the stationary phase method to evaluate the ^-integral in (3.163). Denoting 
the solution of the saddle point equation 



5Seff[A,(t)cl] 

5A 



(3.164) 



by A^'S" and writing A as the sum of ^^f^'^ pi us a fiuctuation, 



A = yl^f^'^ + A, 



(3.165) 



one obtains 



S, 



(3.166) 
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We will now simplify the problem considerably by neglecting the terms of higher than 
second order in A appearing in = — zlndet(P^ + im) |^ Under this condition 

Seff[A,(j)] only contains terms quadratic in A and (3.166) becomes 



Seff[M=Seff[Ait'\M-U[M 

-U[<P,i] = -iln [ Pidet(^^ + im)e-5^^^(^^^)e^/'=''^(-4^^^^^.^). 



(3.167) 
(3.168) 



At this level of approximation we find from (3.162), (3.164) and (3.167) that (j)ci and A 



i) 



are solutions of the equations of motion derived from the action functional 



Siff[A,cP] = S,ff[A,cP] - U[cP]. (3.169) 

The quantum corrections resulting from the path integral over A have produced the ad- 
ditional term —U[(j)]. Under suitable conditions, plays the role of a potential and it 
is worthwile at this point to note some properties of the latter functional. To this end we 

consider an axion field (p = 9 independent of x and denote the integral over the fluctuations 

^ijd^x{-vie]}. 



^jjd^.{-vm} ^ ^^{-um ^ J p^det(^^ + im)e-5^^^(^^^)e^/'^'-(-4^^'^''^-). (3.170) 



V[9] will be called the axionic potential. Changing to Euclidean space using (2.17) and 
pl^ we find 



fd'^xvm ^ I pAdet(^^ + m)e-5^^^''(^^^)e-/'^'-(-4^^-^.^). (3.171) 



Equation (3.171) allows us to draw a certain number of conclusions concerning the form 



of the axionic potential V[9]: 

1. First, we show that e~ '^'^^^^^'^ is real. This can be seen by changing Ao(t,x) 



^o(t,f) = Ao{t,-x), Ai{t,x) A'^{t,x) = -Ai{t,-x), z = 1,2,3 in ( |3.17ip . Un- 
der such a transformation J 9{F A F) changes sign, whereas J F'^ does not. The 
fermion determinant remains unaffected. If 'il}{t,x) is an eigenfunction of for the 
eigenvalue A, then ^'(t, x) = j^ip{t., —x) is an eigenfunction of the D^, for the same 
eigenvalue. 

It is less evident to see whether or not e"-^'^*^^'^^ is positive (and therefore V[9] is 
real). Such is the case for = and we shall admit that V[9] remains positive at 
least for sufficiently small values of 9. 



®Such difficulties may be circumvented by choosing a different approach, see [13]. Since the axion field 
couples to all gauge fields W through a term ~ J 4>{Fw A Fw), one obtains an axionic potential (similar to 
the one which we shall derive) by integrating out additional gauge fields possibly present in the theoretical 
description. 
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2. The last exponential factor in the functional integral (3.171) is obviously real and 
positive. Furthermore the determinant factor det(^^ + m) is also real and positive 

m. 



It was mentioned in section 2.2 that ip^ is a hermitian operator, which there- 
fore has real eigenvalues. Its non-zero eigenvalues are paired in a simple way. From 
{Pf,-f^} = if follows that if ipfil: = Xil), then ipf{-f^^) = -X{j^ij). Hence if A 
is an eigenvalue, —A is also. The determinant det{P^ + m) is then real and positive: 



det(pf + m) = Y[{m - iA) = ]J (m^ + A^) > 0. 



A>0 



Using the fact that e f'^^^^W ig real we find 



(3.172) 



-fd'^xV[9] 



< J DAdet(^i + m)e-/'^'^(-i^'^"'^-) =e-/^'^^M. 



(3.173) 



Therefore V[6] > V[0] and = is the minimum of the axionic potential. 
3. V[6] is an even function of 9. This follows immediately from V*[9] = V[—9] (a direct 



consequence of (3.171) - "*" denotes the complex conjugate) and the point 1 above. 



4. The axionic potential V[9] is a periodic function of 9 and its qualitative shape is 
therefore as shown in figure |3.9[ 



VHMl 




The period of the axionic potential can be obtained from the index theorem [8] 



Using ( |2.28D and ( |2.32D we find 



327r^ 



FAF 



(3.174) 



The ipn are eigenfunctions of the Euclidean Dirac operator corresponding to the 
eigenvalue A„. Since 7^V'n is orthogonal to V'n if A„ ^ 0, only the eigenfunctions 
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corresponding to the eigenvalue A = contribute to the integral in (3.174). This 
leaves 



"3^/^^^ = / ^'^(E4(^)^*(^)-E4(^K(^)) =n+-n_, (3.175) 

\j=l i=l / 

which is the number of positive minus the number of negative chirality zero modes. 
The u[s satisfy P^Ui = 0, jiUi = Ui while P^Vi = and 7^^^ = Vi. This notation is 
consistent with the one introduced in chapter 2. The value of 71,4. — n_ depends on 
the vector potential A, but since it always equals some integer we conclude that the 
period of the axionic potential V[0] is 

Ae = — . (3.176) 



3.4.3 Equations of motion 

The remarks in section 3.4.2 were aimed at explaining the qualitative properties of the 
potential V\9]. Some knowledge of its shape is necessary for finding sensible approximate 
solutions to the equations of motion, which we shall now derive. Recall that by taking into 
account the quantum fluctuations in the A-integral, we have found the following expression 



for the effective action in Minkowski space (see (3.169), (3.168) and (3.158)): 



S'^^^[AA] = W[A]-^ lcPiFAF) + l d''x!^-^Ff^'^F,, + ^{d'^^)id,<P)^-U[<P], 

(3.177) 

with 



W[A] = -ilndet(^^ + im) (3.178) 
-?7[0] = -dn I VAdetiPA + im)e-^^^^^''^^e'^'''<-^^'''^^-^\ (3.179) 

The equations of motion which govern the evolution of the electromagnetic- and the axion 



field are obtained from the effective action (3.177) by calculating 



^^5Z/=o, Ai = 0...3, and ^^^ = 0. (3.180) 

5A^ dcp 

Adding to these the homogeneous Maxwell equations dF = 0, we find the following system 
(^pfiu _ i^fj.upXp^^ denotes the dual of the field strenght tensor) 
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la ~ 5W 
^^^"' = -^29.<f>F''^--JJ (3-181) 

la 

D^ = -^*iFAF)- aU'{4>) (3.182) 
dF = 0. (3.183) 

Explicitly, these equations read: 



^■E = -^V<t>-B-a^-^ (3.184) 



47r2 5Aq 

- - - ■ - - - 5W 

V X B-doE = ^{(PB + V(px E)-ajj (3.185) 
la ^ -> 

D(j)= ^E-B-aU'{(f)) (3.186) 

V • B = (3.187) 

Vx^ + 9oS = 0. (3.188) 
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3.5 Relation between the (4+l)-dimensional and the axion 
field approach 



By comparing the equations of motion written in ( 3.69 )-( 3.73) with those derived in the 



last section, one remarks certain similarities. Thus, before proceeding to the study of 
these systems in the following chapter, we wish to explore the relation between the (4+1)- 
dimensional theory developed in section 3.2 and the system obtained by coupling (3+1)- 
dimensional fermions to an axion field. In particular we will show that the two formulations 
are equivalent if one considers x^-independent fields (more precisely x^-independent vector 
potential A) in the former theory and massless fermions in the latter. In this case the 
component A4 of the (4+l)-dimensional vector potential plays the role of the axion field 



and the thickness of the slab is related to the lenght parameter introduced in (3.150). 



In fact, let us start with the effective action (3.45) derived in (4+1) dimensions 



I "^'"^ {"iL^'''^'^" " ^^''"'^''Af^d^ApdxA-} + rSl^bA], (3.189) 



where 



rSl^l^A] = -Mdet(|)!4U4=i)det(p:4U4=o)} 



(3.190) 



and replace A4 by (p and 84 by 0. Integrating over simply produces a factor of L 
multiplying the first term on the right hand side of (3.190 ). The Maxwell term —j^F^^F^y 
splits up into a contribution —-^F^f^F^p coming from the first four components of the 
vector potential and a kinetic energy term 7^{d^(f)){dfj,(f)) for the field 4> = A^, whereas 
the Chern-Simons 5-form produces a familiar looking coupling between (p and the four 
component gauge field A = (Aq, . . . , A3): 



(3.191) 



^^^^[A, (/)] denotes the effective action for the fields A and = A4 in (3+1) dimensions. 
Since A is x^-independent, the expression for the boundary action rg^[(j4o, . . . , A3)] = 
r^/[(Ao, . . . , A4)|aA] simplifies to 



(3.192) 



an expression, which may be transformed even further. To this end we r e call th e explicit 
form of the operators p\ and P\ in the chiral representation (see (2.42)-(2.43)): 



Pa 



^ 
h 





d 



(3.193) 
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where = a ■ d = —a^do — a-V,^ = a- d = u^Oq — a • V and a subscript A denotes the 
covariant derivative. Inspired by a similar argument in |8] we formally rewrite the product 



of the two determinants on the right hand side of (3.192) as 



det(|)^)det(|)^) = det(-n)det(^^^^) = -det(-n)det(^^). (3.194) 

The prefactor — det(— □) can be ignored, since it merely adds a (diverging) constant to 
the effective action. In this sense we identify 

rJ5J[A] = -iln{det(P^)det(^^)} ~ -ilndet(P^) = W[A] (3.195) 

The effective action obtained from the (4+l)-dimensional theory with x^-independent 
fields therefore reads 



-^F"^F., + i-(a^^)(a,0) 



L 



327r^ 



(p{FAF) + W[A]. (3.196) 



Upon comparision with (3.158) we conclude that this functional is identical to Seff[A,(p] 
obtained by coupling (3+l)-dimensional massless fermions to an axion field. The (4+1)- 
dimensional theory with x'^-independent vector potential and the axion field formulation 



coincide if we identify the lenght parameter / introduced in (3.150) with the thickness L 
of the slab. 

We end this chapter with a remark concerning the physical interpretation of the axion field. 
Our first attempt to devise a (4+l)-dimensional theory was based on an analogy with the 
quantum Hall effect. In that context we defined the quantity {fii — ^ir){x) as being the po- 
tential difference generated by the 4-component of an electric field at the space-time point 
X. The interpretation of ^i^r as chemical potentials of the left- and righthanded fermions 
was suggested by the QH-analogy. In the case of an x^-independent vector potential, 



-LEa 



-LdoA^, 



(3.197) 



see (3.74). Identifying A4 with the axion field (p, we find from (3.197) that the time deriva 



tive of the axion field plays the role of a space-time dependent "difference in chemical 
potentials" between the left- and righthanded fermions. 



Chapter 4 



Cosmic evolution 



The purpose of this report was to present a mechanism which could explain the generation 
of large cosmic magnetic fields in the early universe. In the last chapter we have discussed 
several models from which equations of motion were derived. The present chapter is now 
devoted to the study of these systems of equations. We shall be looking for special solutions 
and try to solve the equations obtained by linearising the system around these special 
solutions. Our hope is of course to find unstable states in the sense that these linearised 
equations predict a growing (electro-)magnetic field. 



We have shown that the axion field theory is equivalent to the (4+l)-dimensional formu- 
lation in the case of an x^-independent vector potential. The equations of motion derived 
from the (4+l)-dimensional theory in turn generalise the system of equations obtained in 



section 3.1. We shall therefore concentrate on the set of equations (3.184) through (3.188) 



Even though we are merely looking for special solutions, a certain number of simplifying 
hypotheses are inevitable. Especially the contribution will be neglected throughout 
this chapter. The argument in favour of this simplification goes as follows: Once the di- 
vergent contribution corresponding to the vacuum polarization graph has been absorbed 
into the Maxwell term through charge renormalization, the remaining contributions are 
of higher than second order in the electromagnetic vector potential A. Thus, if we re- 
strict ourselves to fairly small electromagnetic fields, we may neglect Furthermore, 
if we only consider an axion field that varies slowly in space-time, then we may omit all 
contributions to U[(j)] involving derivatives, dfj_(j), of the axion field. 

The system of equations which we shall consider therefore reads 



V-E 

V xB-doE 
□(/> 



la 

" 4^ 
la 

" 4^ 

V-B = 



i<j)B + V(f)xE) 
E-B-aU'\6] 



V X E + doB = 0. 



(4.1) 
(4.2) 

(4.3) 

(4.4) 
(4.5) 



4.1 Space independent solutions 
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4.1 Space independent solutions 
4.1.1 Solution without axionic potential 

If we neglect the axionic potential U and furthermore assume that (j) is space-independent, 



then equations (4.1)-(4.5) simplify to 



V • ^ = (4.6) 

^ ^ ^ Q. ^ 

\/ xB-doE = --^fiB (4.7) 

dofi = -'^E-B (4.8) 

V • S = (4.9) 
VxE + doB = 0, (4.10) 



where we have set (/) = —jfJ- in order to make it apparent how the system ( 3.75 )-( 3.79) 
obtained in section 3.2 appears as a particular case of the axion field theory. 



A possible special solution to equations (4.6)-(4.10) is F = F^^^ = 



/i = /X 



(0) 



const. 



We may linearise the system of equations around this solution by writing 



B 
E 



ij(0)+,5(l)=g5(l) 
i^(0)+,^(l)=,^(l) 



(4.11) 
(4.12) 
(4.13) 



and ignoring terms quadratic in e: 



V • -E(i) = (4.14) 

V X 5(1) - doE^ = ^(0)^(1) (4.15) 

47r^ 

V • i?(i) = (4.16) 

V X ^(1) + 9oS(i) = (4.17) 

(5o/x«=0). (4.18) 



No particular assumption is made concerning the form of B^^\ E^^"^ and Dropping 
the superscript "(1)" which indicates that we are deali ng wi t h sma ll perturbations and 



introducing the notation fj, 



(0) 



fii — we rediscover in (4.14)-(4.17) the system of equa- 

and charge density 



tions derived in section 3.1 for constant chemical potentials /x^, 
{j^)l3,fl = 0- These equations are linear in the fields E and B with constant coefficients. 
Such a system can be solved by means of Fourier transformation. Using lowercase e and b 
for the (spatially) Fourier transformed fields 
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e(t, k) = j d^xE{t, x)e-'^-^ 
the system to be solved becomes 



h{t,k) = I <rxB{t,x)e 



^\ „—ik-x 



(4.19) 



k-e = 
k-b = 



ik X e + dob = 
ik X b — doe 



a 
'4^ 



{fil - llr)b. 



(4.20) 
(4.21) 



If we choose k = (0, 0, k) in the 3-direction, then it foHows from the above equations that 
only the 2,3-components of e and b can be non-zero. The time evolution of these remaining 
four components is determined by the differential equation 



do 





\ 


/ 










—ik 








\ 
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ik 


4^2 (W -^r) 






62 




bi 









ik 












bi 




\ b2 




V 


—ik 











J 




\ b2 


) 



(4.22) 



The eigenvalues of the above matrix are 

(Eigenvalues)^ = |~^^ =t k^^{fj,i — /ir)| , 



(4.23) 



from which it follows that a real, positive eigenvalue - and hence an exponentially growing 



solution of equation (4.22) - exists for 



a 



(4.24) 



In anticipation of this result we mentioned on several occasion that if in the early universe 
there existed a slight asymmetry in the chemical potentials of left- and right-handed 
fermions (or equivalently a space-independent axion field which was growing at a constant 
rate), this might have lead to the generation of large, cosmic magnetic fields. 



4.1.2 Oscillating axion field and parametric resonance 

Since [/[(/> = 0] for 6 = const is a periodic function of 9 (see section 3.4.2), the space- 
independent solution of the equation □(/> = —aU'[(j)] is either oscillating around (j) = 
or linearly increasing/decreasing with periodic modulations superimposed. In both cases 
the time derivative </) is a periodic function of time, with vanishing mean value for the 
oscillating (j) and strictly positive/negative for the increasing/decreasing fields. 



We therefore consider the equations of motion 



4.1 Space independent solutions 
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V X B-doE 



V -£ = 

la 

V-B = 
V X E + doB = 0, 



(4.25) 

(4.26) 

(4.27) 
(4.28) 



where (f) is space-independent and periodic in time. Again, we will work with the Fourier 
transf orme d vers ions of the above equations and choose k = (0,0, A;). From equations 
(4.27) and (4.28) it then follows that the 3-components of the transformed fields b{x^,k) 
and e{x^, k) vanish. The system of equations for the remaining four components is 



On 





\ 


/ 








Q 


-ik \ 




\ 


62 












ik 


Q 


62 




bi 









ik 








bi 




\ b2 




V 


—ik 








/ 


\ b2 





(4.29) 



where we have put q = — g:^^ for simplicity. The 4x4 matrix in (4.29) - let us call it A 
can be brought to the following block-diagonal form by an appropriate change of basis: 



P-^AP 



( 

-k^- 


V 



1 















1 



,p 



-k'^ + 2kq / 



V 



i{l-k) 


2i 


i(l - k) 


2i \ 


-1 + k 


-2 


1 - k 


2 


-i{l + k) 





i{l + k) 





1 + k 





1 + k 


/ 








(4.30) 



We shall henceforth denote the 2x2 block in the upper left corner of the matrix P ^AP 
by M and the other one by A^. 

At this point we remember that g is a periodic function oi = t. We will now consider a 
particular example and choose q = cos(t), which could be regarded as the solution obtained 
from a parabolic potential approximating T^(</>) near = 0. This approximation is valid 
for small oscillations. 

The problem which we are studying is then equivalent to the Mathieu equation 



+ (a + 26cos(t))y = 0, 



(4.31) 



which written as a system of first order differential equations reads 

do 



M ( ^ 

-a-26cos(t) J \ z 



(4.32) 



Depending on the parameter values, the solution to the Mathieu equation (4.31) is stable 



or unstable. A plot of the stability boundaries as a function of the parameters a and b can 



be found for example in |4J (see also figure 4.1). 
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= 62 + 26 



Figure 4-1: Graphical method for the determination of the instable fc-intervals. The thin lines show 



the stability boundaries of the Mathieu equation (4.31 ) as a function of the parameters a and 6. On 
these boundaries, the solution is periodic. The heavy lines show the intersection of the parabola 
a — {a ~ b'^ + 26) with the instable regions. The corresponding instable fc-intervals are found 
after projection onto the 6-axis, which is identified with k for 6 > and — fc for 6 < 0. 





fc = 0.1 



0.2 



Figure 4.2: The figure on the left shows a stable, but non-periodic solution of the differential 
equation do{y,z) = M{y,z), where M is the 2x2 matrix defined in the text following equation 
(4.30). The initial positions at t = have been chosen {y — 0.7, z = 0) and {y — —0.7, z = 0). 



On the right hand side we plot a solution for a different value of the parameter fc — \k\, which 
illustrates how the corresponding Fourier mode grows by parametric resonance. 



4.2 Space dependent special solutions of finite energy 
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We now compare the 2 x 2-matrix in (4.32) with the 2x2 blocks M and appearing 



in the matrix (4.30), which determines the time evolution of the electromagnetic field. 
The matrix M with q = cos(i) yields the Mathieu equation for the parameters a = k"^ 
and b = k, while the values corresponding to N are a = k'^ and b = —k. The values of 
A; = |/c| for which the Fourier coefficients of the electromagnetic field grow by parametric 



resonance can thus be determined graphically from the plot in figure 4.1 The instable 
/c-intervals are obtained by intersecting the parabola a = 6^ with the unstable regions and 
subsequent projection onto the 6-axis. 

Similarly, for the periodic function q = cos(t) + d (d a constant), which for |d| > 1 could 
be regarded as the solution corresponding to a monotonically increasing/decreasing axion 
field (j){t), one obtains the intervals of instability by intersecting the instable regions in 



figure 4.1 with the parabola a = 6^ + 2db. As before, we have b = k for positive b and 
—b = k for negative b, so the relevant intervals are found by reflecting the (— 6)-axis at the 
origin and taking the union of the contributions from both semi-axes. 

We conclude that in either case the solutions will be unstable for certain intervals of k 
and therefore the magnetic field is growing. The mechanism however is a new one. In 
the example discussed in section 4.1.1, the growth resulted from exponentially increasing 
Fourier coefficients (for small enough values of fc), while this time we found an infinite 
number of intervals, for which the Fourier coefficients grow by parametric resonance. 

It seems reasonable that this general picture remains valid if we consider more complicated 
periodic functions q. In particular, since the instable regions for q = cos(t) were found to 
be intervals and not isolated points, they are stable to small perturbations of the periodic 
function. 



4.2 Space dependent special solutions of finite energy 



Equations (4.1)-(4.5) are Lagrangian equations of motion. They were derived from the 
action functional (3.177) by setting W = 0. The Lagrangian density does not depend on 
time explicitely. Therefore, there exists a conserved energy functional £^ [A, (/>]. The solutions 
discussed in the preceeding section were unrealistic in the sense that they corresponded 
to an axion field (and growing electromagnetic fields) of infinite energy. The instabilities 
in the time evolution of the electromagnetic field are due to a reshuffling of energy from 
axionic to electromagnetic degrees of freedom |13| . 

In this last section we shall argue that the mechanism for the generation of seed magnetic 
fields based on growth by parametric resonance als o work s in systems of finite energy. The 
energy density derived from the action functional (3.177) contains terms in cj)'^, (V0)^, E'^ 



and B^. Hence E, B and d^cf) have to fall off at infinity if the total energy stored in the 
fields is to be finite. This explains the necessity to consider space-dependent solutions. 



4.2.1 Sperically symmetric solutions with magnetic monopoles 

A wealth of special solutions could be derived by allowing for magnetic monopoles in the 
early universe. Writing d^F'^'^ = J^, where Jm = {pm, Jm) is the magnetic current density, 
the system of equations to be solved becomes 
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4. Cosmic evolution 



□c 



Icy 

- — [ct>d^F^'' + d^^F 
la 



327r2 



{FAF)-aU'[(p] 



(4.33) 

(4.34) 
(4.35) 



In vector notation, (4.33) reads 



V-E 

V xB-doE 



la 
~47r2 
la 



{0V ■B + V(f)-B} 



47r2 



{{V xE + doB)(j) + do(l)B + V(/> x E}. 



(4.36) 
(4.37) 



Choosing a sherically symmetric axion field (/'(r), the electric field -E par allel to the mag- 
netic field B and pointing in a radial direction, we can solve equations (4.36) and (4.37) 
by the ansatz 



E 



la 



(4.38) 



Substituting (4.38) into (4.34) yields 



□c 



Pa^ 
levr^ 



B^(j)- aU'[(j)]. 



(4.39) 



Given the form of the axionic potential U, an acceptable special solution for cp must be 
constructed such that 



1. B defined through (4.39) is positive 



2. {(j), B = Bcr, E = Eir} corresponds to a special solution of finite energy. 



However, all of these radial solutions require the presence of more or less strange looking 
distributions of magnetic (and electric) charge, which is the reason why we do not want 
to pursue this idea any further here. 



4.2.2 The sine-Gordon equation and an approximate solution of finite 
energy 

After these purely mathematical considerations, let us recall what we actually intended 
to explain: the growth of seed magnetic fields in a universe with no electromagnetic fields 
initially present. Therefore all the energy is initially stored in the axion field. These physical 
considerations lead us to look for particular solutions corresponding to small (better: 
vanishing) electromagnetic fields and to an axion field which is localized in space. 



We shall try to find spherically symmetric solutions for (j), so that equation (4.3) after 
linearisation around E = B = becomes 



4.2 Space dependent special solutions of finite energy 



63 



9i--^-^)^ = -^U'[^]- (4-40) 
r or or'^ J 

In section 3.4.2 we have shown that U[(j)\ is periodic in (j) (for (f) = 9 independent of x) and 
has a minimum at (f) = 0. In order to work with a concrete example, we choose 

aU'[^] = ^sm{b(P), (4.41) 

where a,b > 0. One idea would be to restrict ourselves to oscillations of small amplitude 
and to approximate ^ sin(60) by a(j), thereby obtaining the equation of motion 

9o - ^ = -^'l^- (4-42) 

Assuming a time dependence of the form e*"^* and setting (/> = ^ we find 

^^ = {a-u')^. (4.43) 
For oj'^ > a there exists a solution which is finite at the origin and falls off like ^ for large r: 



cos 



r 



i{ujt) sin(\/cj2 - ar). (4.44) 



The constant A must be chosen such that Ay/oj'^ — a ^ 1 since otherwise the approxima- 
tion sm(<j)) ~ (/> is not valid. Unfortunately, the total energy associtated with the particular 



solution (4.44 ) is infinite and we have to cut off the function at a certain value of r. The lat- 
ter procedure comes down to introducing a thin double shell of positive/negative magnetic 
(and electric) charge inside of which the axion field is oscillating. 




Figure 4-3: The figure on the left shows the particular solution (4.44) for a/w^ — a = 1 and ujt = tt. 
On the right is a plot of the approximate solution (poHNir ~ 10) for the parameters a = 6 = 1, 
oj = \ and Lot = tt. 
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4. Cosmic evolution 



Another possibility is to keep the potential of the form (4.41), but to neglect the term 



~r'§r Laplacian. This simplification is justified if is appreciable only far from 

the origin. The equation for then becomes the so-called sine- Gordon equation (for a 
spherically symmetric function): 



(9o'- 9,2)0 = -^sin(60), (4.45) 
whose analytic solution can be found in [2lj (see also [3]). It reads 



4 / rjsiniojt) \ \J a — uP 

5DHAr(?",i) = Tarctan — — , r\= , < < V«, (4.46) 

cosh(r/a;r) / a; 



and is translationally invariant (the subscript "DHN" stands for Dashen - Hasslacher - 
Neveu). An (approximate) special solution for (/> which clearly satisfies the requirement of 
finite total energy is therefore 



(\>DHN{r -ro,t) 



(4.47) 



for ro large enough. Of course, (4.47) does not exactly satisfy equation (4.40) for the 



(where the neglected 
4'dhn{'I' — fo) is large). However, it is probably not 



potential (4.41) and needs to be corrected in the vincinity of r 
^ J^-term is large) and r = tq (where 
worth putting much effort into fixing this little deficiency, since our choice (4.41) for the 
functional derivative of the axionic potential was purely arbitrary. 

The purpose of this last section was merely to argue that solutions of finite energy may 
also give rise to growing magnetic fields. We believe that this can be seen from the above 
solution - even though time has not premitted us to attack the difficult problem of sol- 
ving the (linearised) system of partial differential equations with (now) space-dependent 
coefficients. In the neighbourhood of r = ro, where the space-dependence of (p can be 
neglected, we find ourselves in the situation of section 4.1.2, which suggests that in the 
region of space, where the oscillating axion field is appreciable, the (electro-) magnetic field 
will again grow by parametric resonance. 
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Appendix A 



General solution for the continuity 
equation in (3+1) dimensions 



In section 3.1.1 we have mentioned that the the general solution for the continuity equation 
dfj,J^ = in (3+l)-dimensional Minkowski space-time reads 



J° = V-A J = -dtA, 



(A.l) 



where A denotes some vector field and = t. This not so obvious statement will be 
justified below. Assuming that falls off sufficiently rapidly at infinity we can write 



V ■A = r. 



(A.2) 



Equation (A.2) determines the vector field A, whose explicit expression is 



A{x) 



J°{t,y) + B{x), V-B{x) = 0. 



4:7r\x — y\ 



(A.3) 



The first term on th e right hand side of (A.3) is a particular solution of the inhomoge- 



neous equation (A.2), whereas the divergenceless field B{x) is the general solution of the 



homogeneous equation V • ^ = 0. Imposing the condition dtA = — J yields an equation 
for B: 



J{x) = -dtA{x) = -V^ / d'y 



1 



47r|x — y 



V ■J{t,y)-dtB{x), 



(A.4) 



where we have made use of the continuity equation in order to replace OqJ^ by —V • J. 



Equation (A.4) can be solved by means of Fourier transformation. Denoting the spacially 



Fourier transformed fields by lower case j and b, we obtain j{t, k) = T^{k-j{t, k))—dtb{t, k), 

\k\ 

that is 



-^^w k ^ ^ ^ 

dtb{t, k) = j{t, k) - -^{k- j{t, k)) = jtransversal(i, k). 



(A.5) 
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Applying the inverse Fourier transformation to equation (A.5) yields dtB{t,x) 
—Jtia.nsveTsai(t,x), SO the vector field B{x) can be defined as 



B{t,x) 



drJi 



transversal 



(r,x) 



-'transversal(i,^) = J{t,x) - dtV £ / d^y 



1 



Ati\x — y 



■J\t,y). 



(A.6) 
(A.7) 



Since JtransvcrsaiC^, x) is divergenccless we have V • B{x) = 0, as it should. The substitution 
of (A.7) and ( A.6| ) into (A. 3) yields the almost trivial result 



A{t, x) 



dTj{T,x) - V,. / d^y-——j\0,y), 
J 4-K\x-y\ 



(A. 



which proves that any current satisfying the continuity equation can be expressed in the 
form ( A.l ). 



Appendix B 



Gauge invariance and boundary 
currents in (2+1) dimensions 



B.l Gauge invariance 

In (2+1) dimensions, it is possible to write down the explicit expression for the effective 
action. The determinant associated with chiral fermions in (1+1) dimensions has been 
calculated in [20j (for the more general case of a curved space). In this appendix, we 



shall use the corresponding formula for Tg\ (see (B.3) below) and show that the action 
functional obtained thereby is gauge invariant. 

The effective action for the quantum Hall sample discussed in section 3.2.1 is given by 

S,ff[A] = -Scs[A]+rsA[A\9j, (B.l) 

where 



Scs[A] = -^ I AAF (B.2) 



A 



4 



TdA[a = A\QA] = — -r d X ^ aj^{x)a-{x) - a-{x) d y—{x - y)d^a-{y) } . (B.3) 
4 JdA I JdA ^ J 

The following notations have been used in the formula for F^a (recall that the boundary 
of the Hall sample is given by dA = {(xjx^ = 0) U (x|x^ = L)}): 



d+ = do + di d-=do-di □ = d+d- (B.4) 

a+ = Ao + Ai a- = Aq- Ai (B.5) 

and ^{x — y) is the operator defined through the relation 

a,,^{x-y) = 6{x-y). (B.6) 



B.l Gauge invar iance 
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Under a gauge transformation A ^ A + d9 the contribution Scs transforms as 



Scs[A + de]=Scs[A] + ^ [ OF. 

4 JdA 



(B.7) 



The latter result follows from a calculation similar to the one performed in section 3.2.3. 
For simplicity one may suppose that the gauge field 6{x) vanishes on the lower boundary 
of the Hall sample, that is 

e(x)U2=o = 0. (B.8) 
In this case the only contribution to the integral over dA comes from the upper boundary 



X = L and the second term on the right hand side of (B.7) may be written explicitly as 



[ eF = — [ dx^dx^eidoAi - diAo). 

JdA 2 J^2^]^ 



4 JdA 



(B.9) 



The boundary action Tg^ transforms as follows: 



4 JdA ^ 

^ * idld^6) - * - * (dld.e)]. (B.IO) 



a_ - 



The symbol "*" denotes a convolution product. For 9{x) specified in (B.8) this becomes 



T9A[{A + de)\9A]=reA[A\9A] 



dA 



d^x\a+{d-e) + {d+e)a+ + {d+e){d-e) 



\d+e)-{d^e)a^-{d^e){d+e)]. 



(B.ll) 



where we have used = □, formula (B.6) and dx-^{x — y) = —dy^{x — y). The 



expression in between the curly brackets in (B.ll) reduces to 



{...} = -2{Aodie-Aido9). 
After partial integration one therefore obtains 



4 JdA 



From (B.7), (B.9), (B.ll) and (B.13) we find 



(B.12) 



[ d^C{...} = / dx^dx^0{doAi - diAo). (B.13) 



Scs[A + d9] = Scs[A] + ^ [ dx^dx^idoAi - diA^) (B.14) 
Tqa[{A + de)\aA] = TQA[A\dA] + ^ / dx^dx^id^Ai - d^Ao), (B.15) 

2 A2=L 



which proves that the effective action defined in (B.l) is gauge invariant. 
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B. Gauge invariance and boundary currents in (2+1) dimensions 



B.2 Boundary currents 

For an arbitrary A{x), the boundary action F^a yields a boundary current jf^ — 



Using the explicit formula (B.3) for the boundary action of a (2+l)-dimensional quantum 
Hall sample one finds 



J+ = J +J = =a- (B.16 

0O_|_ 

J - /-/ = ^=a+-^*(5ia_) (B.17) 



and the boundary current density is therefore given by 



f = ao-^*{dla-) (B.18) 
/ = -ai + ^*idla.). (B.19) 



Similarly in (4+1) dimensions Tq/^ produces a current density on the (3+1) dimensional 
boundary of the slab. 



